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ABSTRACT 


The Wallops Island Collocation Experiment consisted of an intensive 3-month 
period (April - June^ 1968) of tracking and data collection frOm collocated 
tracking instrumentation at Wallops Islandj using the GEOS -2 satellite as 
target vehicle. • 

The experiment resulted in a wealth of data regarding the ionospheric and 
tropospheric propagation errors, the theoretical and data analysis of which . 
has been documented in some 30 separate reports over the last 6 years. 

This report presents a self-sufficient unifying overview of the entire project 
results with references to the undesdying reports for details. 
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1. DESCRIPTION OF EXPERIMENT 


1.1 Background 

The Wallops Island Collocation Experiment (WICE) (also known as the GEOS -II 
Collocation Experiment) was performed during Aprils May^ and June 1968 as 
a part of the Observation Systems Intercomparison Investigation (OSII), a sub- 
task of the National Geodetic Satellite Program (NGSP). The general objective 
of the OSII program is to improve the accuracy and the estimates of accuracy 
of the geodetic tracking systems through systematic intercomparisons of the 
data afforded by the several tracking systems involved. The WICE Experiment 
in particular was designed to perform certain significant comparisons among 
an extensive array of collocated instrumentation at the Wallops Island Test 
Facility thus effectively avoiding a number of external error sources in the 
comparison due to survey errors, gravitational perturbations, orbital per- 
turbations and timing errors. 

The wide frequency range covered by the assembled instrumentation along 
with the availability of various refraction measuring instruments also afforded 
a unique opportunity for separation and analysis of refraction errors, and a 
considerable effort was devoted to the analysis, development, and testing through 
these means of various proposed algorithms for improved, higher -accuracy 
refraction corrections. 

As the program developed the continuing quest for improved refraction correction 
accuracy led to reanalyses of a number of finer points of thp usual assumptions of 
refraction analysis which are believed to be of some interest for their own sake. 


References 5-13 and 16-21 summarize much of the work done on individual 
tasks under the WICE analysis program. This document is intended to serve 
as an overall guide to the program by expiaining the interrelationships , the 
principal assumptions, -approach and results of the various tasks, referring 
generally to the task reports where possible for details. 
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1.2 Objectives 

The principal objectives of the WICE were; 

• Determination by theoretical and experimental intercoraparisons of the 
practical achievable accuracy limitations of various tropospheric and 
ionospheric correction techniques, 

• Careful examination of the theoretical bases and derivation of improved 
refraction correction techniques as appropriate. 

• Estimation of internal systematic and random error levels of the various 
tracking systems. 

1.3 Organization 

The WICE experiment and its subsequent analyses were performed under &e 
overall directujn of Mr. John Berbert of NASA, GSFC, Experiment Chairman 
and NASA Principal Investigator for the OSII. Other personnel instrumental 
in the final results include 
Experiment 

H.R.Stanley - NASA Wallops - Wallops Project Coordinator and C-Band 
Radar Project Manager 

Dr. H. Piotkin and T. Johnson - NASA Goddard - Laser Representatives 
D. Harris - NASA Goddard - PTH 100 Camera Representative 
P.Kuldeil - Naval Air Systems Command - Tranet Doppler Representative 
R.Vitek, F.Vamura, Dr. F. Rhode, M. Warden - Army Map Service - 
Secor Representatives 

G. Godwin - NASA Wallops - Wallops Project Engineer 
G.Leitar - NASA Wallops - C-Band Radar and BC-4 Camera Data 
Engineer 

T. Savage - NASA Wallops - Timing Engineer 
J.Spurling - NASA Wallops - Meteorological Engineer 
C.Nichols and C. Myers - NASA Goddard - Survey Engineers 
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Data Reduction and Analysis 

H. Parker - RCA - Project Engineer and Technical Information Coordinator 
Dr. A.J.Mallinckrodt - Communications Research Laboratories - 
Refraction Analysis 

Dr. A. Tucker - University of Texas Applied Research Laboratory - 
Tranet Data Analysis 

Reference 5 contains details of the WICE experiment organization and operations. 
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2. DATA DESCRIPTION 


2»i General 

The vehicle for these studies was the GEOS -II satellite whose principal relevant 
characteristics are compiled in Table 2.1. 

The principal data sources utilized in the WIGE analysis include 
NASA Laser Tracker 
BC-4 and PTh 100 Cameras 
C-Band Radars - FPQ-6 and FPS-16 
SECOR Ranging System 

TRANBT System (Including the 3 -frequency TRANET) 

Ground Meteorological Instrumentation 
Balloon Sondes 

Bottomside Ionospheric Sounding 

(Wallops, Grand Bahama and Ottowa) 

Topside Ionospheric Soundings 
(Allouerte) 
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TABLE 2.1 


GEOS n CHARACTERISTICS 


Orbit: Altitude 

Pd 

l£2i km 

Eccentricity 

Ad 

.03 

Indination 

Ad 

105° 

R. A.Asc Node 

Ad 

275° 

Avg. Perigee 


260.3° 

Mean Anomaly 


83.6° 


As of Epoch C/25/1968 


Satellite Instrumentation 

C-Band Transponders (2) 

C-Band Passive Reflector (Retrodirective Array) 

U.S.Navy Doppler Beacon 

Comer-Cube Optical Reflectors 

Optical Beacon 

NASA Minitrack Beacon 

NASA Range and Range -Rate Beacon 

U.S. Array SECOR Beacon 
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2,2 Laser 

The laser tracking system, which was moved to Wallops Island for the duration 
of this test, was built and operated by the Optical Systems Branch (OSB) of 
Goddard. This system uses an intense, highly collimated, short -duration beam 
of light for illuminating the spacecraft being tracked. At the spacecraft, the 
beam is reflected back towards the ground station by an array of cube comer 
reflectors. The returning light is detected photoelectrically, and its transit 
time is measured to yield the range data. The actual laser transmitter is 
mounted on a radar pedestal along with a Cassegranian telescope used for 
receiving the reflected laser beam. When the laser system is tracking, the 
transmitter is flashed at 1 pulse/sec. Each transmitter pulse starts a time- 
interval -measuring unit necessary for range measurement. During the pass, 
the mount, equipped witli digital encoders, is directed toward the expected 
position of the spacecraft by a programmer fed with punched paper tape. By 
using a telescope, the operator can see the spacecraft and make corrections 
to keep it witliin the illuminating beam, which is only about 1.2 milliradians 
wide. Along with a range measurement, both the azimuth and elevation of the 
spacecraft are recorded from the position mount. At the time of the.se tests, 
in 1968, the laser tracking system was probably unbiased to 0.15 meter in range, 
with a random noise component of about 1.2 meters, and could produce range 
rates through an orbital fit to the range data which were good to about 1 cm/sec. 
Now, in 1975, the laser bias and random noise have both been reduced by an 
order of magnitude. These estimates include aU known error sources except 
the scaling error of 1 part in 10^ due to the uncertainty in the velocity of light, 
which affects ail systems . Mount angular measurements are recorded but are 
considered as a secondary measurement since they are dependent on manual 
tracking. 


2.2.1 Laser Data Preprocessing 

The OSB personnel were responsible for laser data preprocessing. The pre- 
processor program accomplishes the following functions: 

- Converts the recorded time of observation to the time when the laser pulse 
was at the spacecraft. 
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• Computes the range to the satellite from the round-trip time interval 
values and calibration values. 

• Corrects the measured elevation an^e for refraction. 

• Corrects the computed range for refraction. 

• Edits the data based on a five -sigma rejection criterion. 

• Reformats the acceptable data points into the required Geodetic Satellite 
Data Service (GSDS) fosnnat and outputs the data on a magnetic tape, with 
a density of 1 observation per second. 

Preprocessing details are contained in Reference 14 and Reference 15. The authors 
received the WICE laser data from GSES and conducted this intercomparison study 
with no additional preprocessing. 


222.2 Laser System Calibration 

For angle calibrations, a special boresight feature is incorporated in the collimating 
optics for the transmitted beam, which allows the laser transmitter to be aligned 
parallel to the opt -mechanical axis of the tracking pedestal. Boresighting is 
accomplished fay firing the laser through a separatdLy attached focusing lens onto 
a piece of aluminum foil in its focal plane. The reflex viewer, which forms the 
boresight function on the collimating optics, is then inserted in the optical path, 
and its cross-hairs are adjusted to coincide with the hole formed in the foil by the 
focused laser beam. With the focusing lens removed, the reflex viewer is directed 
along the laser beam and can be used to bring the laser optical axis parallel to the 
other optical systems on the tracking pedestal. 

For range calibrations, the total delay in signal due to td.es cope optical path length 
and delay through the photomultiplier tube is measured over a geodimeter calibrated 

range (3274.98 meters at Wallops) before and after each pass. See Reference 14. 

1 
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2.2.3 Laser System Timing 

The laser data control unit generates all the control signals for operation of the 
laser and receiver systems. In addition, the unit maintains system time with respect . 
to an external time source such as V/WV or, as in the case of WICE,to the * 

Wallops station master clock. This is accomplished by setting the laser control 
unit. A l-MHz oscillator, acting as a secondary time standard, is counted to one 
pulse per second through phase shift and delay circuits for synchronization with the 
external timing standard. At WIGE the laser 1-pps signal is s 3 mchroniz;ed to about 
+ 0.05 millisecond prior tc each pass with the master clock cable signals adjusted 
for a cable delay bias and the current delay between the Ce standard and the TODG. 

The l”pps signal is then used as the on-time generated pulse throughout the entire 
data control unit and operates a binary coded decimal (BCD) time code generator • 
whose output is displayed visually as well as recorded through the data select gates 
for correction witli measured range. The rotation prism Q-switch cannot maintain 
exact synchronism with the on-time pulse, and therefore, the laser may fire from 
8.5 to 11 milliseconds after the command time. An uncertainty of this magnitude 
in the time of observation is not compatible with the accuracy requirements, so a delay 
time interval counter was incorporated in the data control unit to accurately 
measure the time of firing with respect to on-time. This counter is started by the 
on-time and stopped by a signal from the laser beam sample unit, giving the 
absolute time at which the laser fires to within 100 microseconds. The output of 
the d^y counter is stored and transferred to die data select gates in parallel with 
the range-time-interval measurement for recording. This value is used in the 
preprocesser to correct the data time tag. 


2.2.4 Laser System Tracking Constraints 

The GSFC laser located at Wallops had the following tracking constraints: 

• Nighttime at the station (sun 10° below the horizon). 

• Satellite maximum elevation angle above 30°. 

. Satellite sunlit or lamp flashing for visual acquisition. 

• One safety operator at laser station to make visual observation for 
low -flying a ir craft . 

• An operational surveillance radar to verify that there were no aircraft 
within a 14 -nautical -mile range. 
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2.3 SJBCOR 

The SECOR system, developed fay the CUBIC Corporation provides a highly accurate 
modulation ranging multi -lateration system with inherent ionospheric range error 
estimation and correction. SECOR was operated for WICE and the data preprocessed 
by the Army Map Service (AMS). 

The system utilizes carrier frequencies of 420. 9 MHz up and 224. 5 and 449.0 MHz 
down. The difference in measured range on the two (2:1 related) down frequencies 
provides the ionospheric error estimate. Range is measured by a series of sine 
wave modulation tones starting with 585.533 KHz and ranging down to 20 Hz for 
ambiguity resolution. System resolution is 0.25 meter. The geometry is pure 
ranging multilateration, each ground station interrogating the transponder 
sequentially and utilizing only ti:e downcoming response to its own interrogation. 

f 

2.3.1 SECOR Data Preprocessing 

AMS personnel were responsible for SECOR data preprocessing. The prepro- 
cessor program accomplishes the following functions: 

• Computes the time of observation, which is defined as the time when the 
pulse was at the spacecraft. 

• Makes ambiguity corrections to the edited range measurements. 

• Applies calibration values to the edited range measurements . 

• Applies tropospheric refraction to the range measurements. 

• Uses the difference between ranges measured on the low - and high-frequency 
carriers from the spacecraft to compute a correction for retardation due to 
the ionosphere. If 

= measured range at f^^ = 224.5 MHz 
= measured range at f„ = 449 . 0 MHz 
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then^ since the ionospheric range errors are to first order proportional 
to 

= R + Kf~^ 

M2 = R + Kf2^ 


which can be solved for 


R 


and 


K 


or 


AEj 


^Mi -4 m2 



ionospheric error on f 2 


M 2 - R 


K£, 


-2 


M^ -M 


T 

I 


- 1 


2 


The range is corrected for this value, and the ionospheric correction value 
is included in the output. 

• Reformats the data into the required GSDS format and outputs the data on a 
magnetic tape, with a density of 1 observation per 4 seconds. 


■J 


2.3.2 SBCOR System Calibration 

In the calibrate mode, the calibration oscillator generates 196.4 MHz. This is 
fed to a mixer mounted above the vertical axis of the WICB station dualreflector 
antenna system and between the up-link and down-link reflectors. A 420. 9 MHz 
ground station up -link carrier frequency is radiated, from the up -link antenna 
to the mixer, to produce 224.5 MHz and its second harmonic, 449.0 MHz, which 
are the spacecraft down-link carrier frequencies . These are reradiated to the 
down -link antenna, providing a closed-loop method of determining approximate 
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I i 

; i 

zero-set of the range servos prior to each pass. All components in the ground ' 

station are inside the calibration loop. j j 

A refinement of the zero -set is made by air link calibration prior to each tracking 
pass and immediately after each tracking pass . This utilizes a mixer as above 
and a discone antenna which is 28 meters from the SEGOR su 2 rvey reference mark^ 
and is fed through a cable connected to the ground station. 

The difference between the 2 8 -meter surveyed range to the discone and the 
measured range is recorded on calibration sheets, for both the high- and low- 
frequency channels, for doth precalibration and postcalibration measurements. 

These calibration numbers are used in the preprocessor to correct the range data. 


2.3.3 SEGOR System Timing 

The WICE SEGOR station has a rubidium clock. The rubidium clock was used 
to operate the time code generators which record UTG time on the magnetic tape 
with a resolution of 1 millisecond each time the digital seorvos record the range 
on the tape. 

The Wallops Island range time was derived from an HP 5060A cesium clock, set 
to UTG OSIAVOBS). This clock was periodically transported to the SEGOR site 
in order to check the rubidium clock. The offset between the rubidium clock and 
UTG ONfAVOBS), as recorded on the data logs, was always between 5 and IS micro- 
seconds during the WICE operation. 


2.4 TRANET 


The TRANET system transmits a set of accurately determined frequencies from 
the spacecraft and provides one-way doppler measurements to the ground station. 
TRANET utilizes three harmonically related transmitted frequencies of approxi- 
mately 162, 324, and 972 MHz to provide ionospheric corrections an<^as well. 
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corrected data (doppler frequency) from any two of the diree. TRANET was 
operated during WICB by University of Texas personnel under the direction 
of the Johns Hopkins Applied Physics Laboratory (APL) and the Naval Air 
Systems Command. The quantity measured by a TRANET tracking system is 
the doppler frequency as a function of time. 


2.4.1 TRANET Data Preprocessing 


The Naval Weapons Laboratory (NWL) personnel were responsible for the 
TRANET data preprocessing after collection of the data by APL. The TRANET 
data underwent the following preprocessing; (Ref. 15) 


• First-order ionospheric refraction error estimates and corrections are 
made by analog techniques, using equipment at the tracking station. This 
correction is based on the assumption that the ionospheric doppler error 
varies inversely ivith the carrier frequency so that die corrected doppler 
shift and the doppler error referred to fj^ are given by 

f f f “ f^f 
^1^2^02 1 D 2 

^ ^2 1 



• f “ f f f 
•2 U^2^D2 . 

^2 h 


• The time of observation is computed. This is defined as the observed time, 
at the station, of the midpoint of the doppler integration interval. The 
calibration value (offset of the TRANET station clock from the Wallops 
Island cesium clock) is used to correct the observation time to UTC (NAVOBS). 

• The observation frequency is corrected for the- error in the station frequency 
standard determined from VLF comparisons. 

• A spacecraft reference frequency (base frequency) is computed for the pass. 

• The data are edited based on a 2.5-sigma rejection criterion. 
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• The remaining observations are aggregated in groups of eight, covering 
a 32-second interval, A smoothed frequency value is calculated by fitting 
a straight line to the residuals in the 32 -second interval and evaluating 
the fit at the central time of the interval. The residual corresponding to 
the fit at the central time is then added to the computed frequency for that 
time. These data are run through a reformat pregram which arranges the 
filtered data into the format required by the GSDS. All smoothed frequency 
values and the base frequency for each pass are scaled to 108 megacycles 
by multiplying by ■■ , where fe is the nominal equivalent fre- 

quency obtained from a table and approximates the frequency out of the 
station refraction corrector unit. 

The TRANET are the only type of data which undergo mathematical smoothing 
in these intercomparisons. 


The WICE TRANET data was further processed at Goddard through the following 
additional preprocessing steps prior to intercomparisons with the other systems. 
- Conversion of the recorded time of the observation from observed time 
at the station to the time the signal was at tlie satellite by subtracting 
one -half of the round trip time. 

• Conversion to range rate values in meters/sec by the following 

® 

where Fg = Base frequency received from NWL 

= Smoothed measured frequency received from NWL 

c = Velocity of light 
“ 2.997925 X 10® m/s 


• Correction for tropospheric refraction, using the formula 
AR^ = 


8432,336 N EGosE 
s 

X" 


(SinE + .026) 


i 


■] 

1 
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A = correction (cm/sec) to add to R 

Ng = ground refractivity, (U - X) 

E = elevation angle 
B = elevation angle rate (radians /sec) 


2.4.2 TRANET System Calibration 

The station frequency error which appears in the doppler data header in the 
teletype to APL is the departure of the frequency of tlie station standard as 
determined from a known (VLF) refesrence frequency . This known correction 
is applied^ in the NWL preprocessing program, to the frequency measure- 
ments. 

A nominal value of ttie satellite oscillator frequency is associated with each 
spacerraft but is modified for each pass as follows. First, NWL computes 
0-G*s by comparing the VLF corrected doppler frequency measurements with 
the doppler frequencies predicted from a reference orbit. The reference orbit 
is determined with previous doppler data from the entire TRA?'FT network. 

The 0-C's aie then used to compute an estimated frequency bia^ for each pass. 
The spacecraft oscillator nominal frequency corrected for this bias is called 
the base frequency and is included, as an additional number, with the frequency 
measurements submitted to the GSDS for each pass. 

Since the deteomiination of the base frequency involves the entire TRANET 
network, the WICE TRANET data are influenced by this network. Data sub- 
mitted from the other WICE systems are not influenced by any other stations. 


2.4.3 TRANET System Timing 

The station clock error accompanying the doppler data in the teletype to APL 
combines the station clock offset from the received Wallops Island pulse, the 
cable delay, and offset of the Wallops Island working clock (TODG) from the 
cesium dock. The doppler data submitted to GSDS are referenced to the 
Wallops Island cesium dock which is set to UTC (NAVOBS). 
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2.4,4 Three -Frequency Data 

In addition to the above described "standard" TRANET data there were made 
available thx'ough the University of Texas, separately received and recorded 
three -frequency "Geoceiver" data. From this data it is possible to derive 
explicitly both the ionospheric corrected data and the ionospheric correction, 
per se, by the equations given above. 


2.5 C-Band 

Data were received from two G-band radars operated by Wallops Island, the 
AN/FPQ-6 and the AN/FPS-16. Both are pulsed radars capable of non- 
ambiguous range measurements of up to 32,000 nautical miles, and each provides 
azimuth and elevation angle measurements to the target. The FPQ-6 radar 
can also measure range rate (R) if used with a coherent transponder or if 
the reflected signal from the spacecraft is strong enough for skin tracking. 

A passive retro -directive Van Atta array on GEOS-2 makes it possible for 
the FPQ-6 radar to skin -track this spacecraft. 

Two C-band beacons were installed in GEOS-2. Beacon #1 has a 0.7-|jLSec 
fixed nominal delay, and beacon #2 has a S-|isec fixed nominal delay. 

The FPQ-6 radar’s subsystems may be functionally grouped under signal 
detection (transmitter, antenna, and receiver), target acquisition, target 
tracking (range and angle servos), data processing, and system control. An 
ultrastable frequency-synthesizer-multiplier chain , power amplifier, and 
hard -tube modulator form the C-band transmitter. The antenna comprises a 
solid -surface 29 -foot parabolic reflector illuminated by a monopulse, polari- 
zation-diversity cassegrainian feed. This structure is supported by a 2 -axis 
(azimuth -elevation) pedestal featuring a low -friction hydrostatic azimuth 
bearing, anti-baddash drive gearing, and precision single-space 20-bit - 
angle-shaft encoding subsystem. The angle, or antenna -positioning, subsystems 
are high torque -to -inertia electrohydraulic servo loops. Tracking signals are 
supplied to the antenna-positioning and ranging servos by a low-noise, broad- 
band, 3-channel receiver subsystem. An all -electronic digital ranging sub- 
system affords unambiguous range coverage to 32,000 nautical miles at high- 
pulse-repetition rates, with a granularity of 2 yards. The data system con- 
tains a 4096 -word coincident -core, bus -organized, stored-program, militarized 
computer (RCA, FC-4i0i). 
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The FES -16 radar is very simRar to the FPQ-6 except tliat it has a 17 -bit 
an^e encoder and a 12 -foot parabolic reflector and does not have a computer 
and skin track capability. 


2.5.1 C -Band Data Preprocessing 

The following preliniinary preprocessing is done by Wallops Island: 

The on-site RCA 4101 computer program for the AN/FEQ-6 was used to apply 
the static corrections (pedestal mislevel, droop, nonorthogonality, encoder bias, 
encoder nonlinearity, and skew) to the raw FPQ-6 data, but not to the FPS-16 data. 

Dynamic lag corrections calculated by the 4101 program are recorded, but are 
not applied to the data. The 4101 FIS -16 raw data tapes were processed through 
the Wallops preprocessing program which applies a time tag correction to the 
data,. converts the data from radar bits to range in feet and azimuth and elevation 
in decimal degrees, and reformats the data from 4101 format to the standard 
GEOS-B radar data format, sometimes called the modified Calsat format. 

WICE - C-bahd data in the modiEied Calsat format were sent to the Principal 
Investigator, GEOS OSII. These data were then additionally preprocessed by 
the WIGE-G-band preprocessor program at Goddard, which does the following: 

• Computes the time of observation, which is defined as the time when the 
piilse was at the spacecraft. 

• Applies tropospheric refraction corrections to both the range and angle 

measurements . f 

• Reformats the data to a format compatible with the GEOS data adjustment 
program (GDAP). (Editing is done by hand after residuals are obtained with 
GDAP against the laser reference orbit.) 

• Selects every Nth point (one per second). 

• Applies range bias correc»._ons derived from the appropriate nominal 
beacon delay and from the pre- and post-pass range target measurements. 

% 
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2 , 5 o 2 C -Band System Galibratioa 

For pre- and postmission calibration, data tape recorders are run for approxi- 
mately 10 seconds (at 10 samples/second, this gives approximately 100 samples), 
recording each of the following: 

• Selected AG C values. 

• Boresight tower (BST) normal - Antenna electrically locked to the BST in 
azimuth and elevation. 

• Boresight tower plunged - Same set-up as for BST normal, except the antenna 
is in the plunged mode. 

• Range target, skin gate - (Ef transponder track is planned, die proper delay 
compensation should be set into range system prior to this step.) Lock on 
range target in skin gate. Range displays and recordings should read the 
surveyed range to the range target. 

• Range target, beacon gate - Range displays and recordings should read the 
surveyed range to the range target minus die proper delay compensation. 


2, 5. 3 C-Band System Timing 

The received pulse from the Wallops master clock is used by the G-band systems 
to time-tag the range, azimuth, and elevation data. The circuitry and cable 
transmission delays bias the time tag by +5.90 milliseconds for the FPQ -6 
and +1.05 milliseconds for the FIS-16. These known biases are accounted 
for in the Wallops program by adding them to the recorded data time tags to 
give UTC time. 

The delay between the Wallops Ce Standard and TODG (200 + 100 lasec) has not 
been included thus far in the G-band data time tag corrections. However, these 
variations are measured and recorded daily at the master site by direct compari- 
son of the TODG oscillator and the cesium beam standard, so the proper cor- 
rection to the time tag can be applied at some future date. 
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2,6 Cameras 

Data -were taken on GEOS -2 flashes from four BC-4 cameras operated by 
Wallops Island and one PTH-100 camera operated by Goddard, These data 
were utilized in addition to the laser data as a primary reference in the 
derivation of GEOS reference orbits. 


2.7 Meteorological Instrumentation 

Ground measurements consisting of Temperature, Pressure, and Relative 
Humidity were recorded from measurements at the FFQ-6 radar site for each 
pass. Estimated accuracy of these measurements is 

• Pressure + .01 " Hg 

• Temperature +1.0°F 

• Relative Humidity + 1% to + 3% 

Radiosondes were released for P,T, and RH profiles within 10 minutes of the 
start of each pass. Altogether some 93 such profiles and associated ground 
level measurements are available. These have been reduced to refractivity, 
ray traced, and utilized as the basis for various comparative studies discussed 
in section 4. 


2.8 Ionospheric Soundings 

Routine periodic (generally 15 minute) bottoms ide ionospheric soundings were 
available throughout the WICE experiment from stations at Wallops Island, 
Ottowa , and Grand Bahama. 

In addition, near coincident Alouette 1 and 2 topside soundings were avaUabie 
for some passes at randomly related time and positional (Point-of-Closest- 
Approach) differences relative to the GEOS -2 passes. In order to select reason 
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ably "coincident" passes a somewhat arbitrary meas nre of correlation was 
defined in terms of the expected time and positional correlation scales and 
discrepancies as follows: 

p = exp(-AT/2) exp(“AR/3000) 

where AT = time difference of POCA, hours 
AR = distance between POGA’s, miles 

advanced or retarded'to common time at sun line rate. 
Differences are between GEOS -2 and either Alouette 1 
orAlouette2. 

17 passes were thus identified having p > 0.4 and these were taken as the 
standard ionospheric test cases for future analyses. These are identified 
in Table 2-2. 


r 

r 
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3. RAY TRACING 


REBK Program . ^ 

Since it is tised as a standard in many of the refraction comparisons to be used 
in ensuing sections we describe in this section the basic ray-trace program, 
REEK (Ref. TrimOle 1970), used in these studies and later some investigations 
of fine points with regard to its use. 

The problem addressed is that of determining the phase path, its bending, and 
phase path length in a spherically symmetrical refractive medium with arbi- 
trarily specified refractive index vs height profile (assumed isotropic). 


REEK solves the differential equations for the phase path (wavefront normals) 
in the form: 


de 

Hi- 

_ COScp 
ir+h 


(S.1-1) 

dh 

Hi- 

= sincp 


(3.1-2) 

de 

HH 

1 % 
■ 5“ "ar 

P 

(3,1-3) 

dp 

Hi- 

= lL 

P 


(3.1-4) 

cp 

= 0 + 6 



where (see Figure 

3-1) 



s = 

geometric path length 



6 = 

earth central angle subtended by path 



P = 

effective radio phase path length 



h = 

height 



6 = 

elevation angle of ray relative to horizontal at takeoff point 


“p ^ 

phase refractive index 



cp = 

local elevation angle of ray 


These equations are solved by numerical methods, subject to the following 


boundary conditions: 



0(0) 

= 0) 



h(0) 

= 0 ( 

starting at ground receive r 


p(0) 

= of. 
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9(s.) = 0 1 

y target coordinates 
h(Sj) = U 

The last two are equivaleat to saying that for some initially unspecified^, 
the ray should pass through the prescribed target coordinates, REBK inputs 
are in terms of either true or apparent elevation angle or range. In the case 
of apparent inputs the solution is straightforward starting the ray with initial 
conditions corresponding to apparent elevation angle and integrating until 
apparent range is reached. In the case of true inputs ^ the program uses a fast 
converging iterative method to find the initial ray angle corresponding to the 
specified end point coordinates. ' * 

To achieve maximum numerical accuracy, the equation for p , given in eq, 3. 1 - 
is decomposed by defining 

Bp = p-R 

where Sp = total phase path error 
R = true range 

Bp = 

®B ®R 

Where the first term is defined as 

Bg — s - R = phase path bending error 

= p -s = phase path retardation error 

end both Sg and are small quantities and can therefore be numerically 
integrated more accurately. Note that the solutions for the bending and retarda- 
tion terms could be written, once the path (P) is found, by solution of eqs .0.1-1) 
(3.1-3) 



(3.1-5) 


(3.1-6) 

(3.1-7) 


(3.1-8) 

(3.1-9) 
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One minor modification was incorporated in the REEK stibroutine for the present 
program. In the original version, refractlvity in or abpve the ionosphere was 
extrapolated linearly which in some cases resulted in refractivity going through 
zero to Hie opposite sign. In order to improve the accuracy and eliminate any 
such problem, the routine was changed to use exponential interpolation or 
extrapolation above the top of the ionosphere (Ref. 21). 


3,2 Straight Path Approximation 

For some purposes it is reasonably accurate to ignore bending and approximate 
the ray path by the geometrical straight line from transmitter to receiver. The 
computations under this approximation are greatly simplified. In order to gain 
some understanding of the limitations of this approximation and for other uses to 
be developed later it is useful to consider a simplified ionospheric model for 
which the exact and straight line approximate solutions can be written analytically. 

For this purpose the ionosphere is modelled as a simple planar uniform slab as 
shown in Figure 3.271 where 

h “ height of sateEite 
T - thickness of ionosphere slab 
N = refractivity (phase) of ionospheric slab = n-1 
E = 90° -cp = elevation angle geometric line-of-sight. 

The exact analytic solution for bending and range errors can be developed for 
tfam model (Ref. 6). If we let 

denote geometrical path length of some path x 

p denote refractivity integral along some path x 

A 

S subscript^ denotes straight path 
B subscript^ denotes bent path 
p,g subscript, denote phase or group respectively. 

Then we can define four types of "range" of interest 
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R = J ds = true geometrical range 
S ' 


p g » J* ”-p ' radio range using the £hase refractivity 

P S integral along the straight path . 

R_ = J* ds = Geometrical length of the bent £hase path 

P Bn 




p „ = J n ds = radio -range using the phase refractivity 

P P Bn ^ 


integral along the bent £hase path. This is 
•what the measuring system actually senses. 

p g is, of course, the actual measured radio range while p g is an often 

P P P 

used approximation. Rg is equivalent to s, the variable of integration in 

P 

the REEK formulation 


Phase Path 


From Snell*s law and the geometry, approximate solutions for the bending angles 
defined in Figure 3.2-1 may be developed as power series in N; 


P = -NtanS +N tanS 




+ 0 , 


a = -^Ntan® ^tan$ 


1+ |tan§ - J-{i+2tan^§) 


+ 0 , 


where O. = 0(N^) denotes neglected terms of order and higher. 
J ^ 


Then for the various path integrals as defined 


(3.2-1) 


(3.2-2) 

(3.2-3) 
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R = hsec$ 


(3.2-4) 

p S - R 

P 

+ NpTsec® 

(3.2-5) 

11 

IS[?T T o 

+ — ■* U") seel +Og 

(3.2-6) 

Pg = R 

P P 

N?t 2 

+ NpT seel ^ — (1 - ^) tan I seel + Og 

(3.2-7) 

It is interesting to 

note that the error due to ignoring bending is 


(P B “ P = 

P P P 

-(Rg -R) 
P 



■ 9 

— — (1 - tan s sec$ 

(3.2-8) 


i. e . , the bending reduces the measured radio path length_, p g j relative to 

P P 

the straight path integral by just Ihe same amount that it increases the geo- 
metrical patli lengthj and this reduction or increase is strictly a second order, 
term. The dependent term may legitimately be called a bending term 
since it occurs only in connection with the bent patlis and reduces to zero for a 
case of normal incidence (? = 0) where there is no bending. 

Group Effects 

It is possible to carry one interesting step further and examine group effects 
in the same model for, quite generally, the group measured range is related 
to the phase measured range by 

For the purposes of this example we will assume that the ionospheric refractivity 
-2 

follows a perfect f dispersion law. Accordingly 

Np = -Kf-2 

and from (3.2-7), recognizing that p = pg 

P p p 


(3.2-10) 
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3ivt2 _ /t T\ * — 2 , 


Pg - R“NpTsec5 + jNpT(l - tan fisecf+Og 


Note, interestingly, . the same resiilt is obtained if we were to integrate 
the group refractivity along the phase bent path, as is commonly done, for 
example in REEK, i.e., 

P B “ J “s'*® 


g P B, 


= J(l+Ng)ds 


Br 


= J(l-N)ds 

= R“NpTsec5 +^NpT(l -^)tan^$ sec5 +Og 


which is identical to (3.2-12). 


Similarly we can derive 

r, - d 
" df 


g 


(fp_s) 


= R-*NpTsec8 


and thus the difference between the straight and bent path computations is 



This difference is so small that it can often be ignored for practical purposes. 
In terms of the ratio of errors computed in the two ways 



(3.2-12) 


(3.2-13) 


(3.2-14) 


(3.2-15) 


(3.2-16) 
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To illustrate the error in the straight path assumption with a numerical 
example for the ionosphere, consider a case where 


so 


T = 293. 57 km 
h = 1333 km 

f = 434 MHz (SECOR effective two-way frequency) 
f^ = 5.653 MHz = vertical incidence critical frequency 



Then as a function of 5, we have from equation (3.2-16) 


Elevation Angle 
at Ionosphere 
B., Degrees 

90 

45 


AR 3 /AR g 

1 ^ (vertical ray - no bending) 

1.000099 


20 


1.000749 


There is no need to consider E. less than 20°, since the elevation angles 
in the real spherical ionosphere cannot be much less. 


As a further check of the straight path assumption in a less idealized case a 

numerical comparison of the straight and bent path integrals for a spherical 

earth was carried out for a Chapman ionosphere with a maximum refractivity 

“A 

(N = -84.84 X 10 ) and the satellite at 1333 km. The straight path integral 
was carried out by a specially developed straight line raytrace program using a 
simple trapezoidal integrator and the bent path integral was carried out by the 
REEK raytrace program, modified for group range errors. The resulting dif- 
ferences between the outputs of these two programs are presented in Table 1 . 
The difference, at most, is 0.11 meter out of 73 meters total refraction, or 
0.15%, and is in reasonable agreement with the simple theory of the difference, 
equations (3.2-15,16), at low elevation angles, but is dominated at the higher 
elevation angles by a bias of about 0.04 m, presumably due to difference in the 
REEK and straight line integration formulas. These results are taken as 
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coofirming the approximate equations (3.2-15 and (3.2-16). 


Table 3.2-1. 


Difference Between Straight Line Raytrace (AR g) and 

REEK Raytrace (AR g) ® 

g 


Elevation 

(degrees) 

(E-90-4) 

ARg 

g 

(meters) 

AR g - AR g 
g g 

(meters) 



Numerical 

Theoretical 
(eq. 3,2-15) 

O.T 

73.336 

-0.1140 

- 0.0621 

1 

73.245 

- 0.1122 

- 0.0618 

10 

65.812 

- 0.0715 

-.0.0432 

20 

52.960 

- 0.0424 

-0.0199 

40 

35.911 

- 0.0355 

- 0.0039 

60 

28.311 

- 0.0395 

- 0.0008 

80 

25.360 

- 0.0428 

- 0.000069 

89 

25.027 

- 0.0434 

- 0.0000006 


The above is based on a spherical earth and a Chapman profile where 

N = N^^^exp(l -z -e 
max 


h“h 


= -84.84 x10", 
max , 

h„^^ = 375 km 
max 

Hg = 108 km 

h = 1333 km 
sat 

The figures in the theoretical column are based on equation (3.2-15) corrected 
for angle of incidence at 375 km^ and using 


T 


108 e 


293.57 km. 
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For the troposphere with 

T = 7 km 

h„ = 1333 km 
s 

N = 350 x10 

we have from equation (3.2-15) (using the positive sign for Np 
troposphere) f ■ ~ ! !' 


,-6 


iilevation Angle 
degrees 


AR „ 

g® 

(meters) 


AR g-AR g 

g g 
(meters) 
(eq!n. 3,2-15) 


in the 


AR g 

—-S-. 
AR <3 
g® 


2 

100.26 

30.06 

, .7001 

3 

55.71 

. 8.90 

1 .8402 

4 

38. 87 

3.75 

.9035 

5 

30.02 

1.91 

I .9361 

10 

14.34 

.23 

i .9834 

20 

"7.19 ■ ! 

.03 

' .9960 


Thxis for accuracy of the order of 2% or .2 m which is near the limit of correction 
accuracy for the troposphere the straight line assumption is limited to elevation 
an^es of 10° or greater. This is a severe limitation on the applicability of any • 
approach ignoring bending in the troposphere. 


3,3 Group vs Phase Effects 


Basic Relations ’ ' 

' ' ■ — r 

0 

In the ionosphere j the index of refraction is a function of frequency, i.e. , 
dispersive and modulation or group eSects travel at a different velocity and 
are subject to different overall delays than phase effects. 


In general, between any two points in a linear network, die group and phase 
delays are related by the classical relation 

'^g = 

which may be taken as defining the group delay. 


(3.3-1) 
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Since inferred range (group or phase) is proportional to t (p = Tc), then, 
irrespective of the inhomogeneity or path taken in the intervening medium, the 
same relation holds between ^oup and phase range 

d 


'’g = 


df 


(fpn) 


(3.3-2) 


In principal, this provides a precise basic means of determining group range 
error from any accurate ray tracing program by numerical differentiation of 
the terminal, i.e. , end-to-end phase range error results as a function of frequency. 
In practice, however, this involves several times as much computation as a 
single ray trace and is subject to the numerical problems characteristic of 
numerical differentiation so a more direct means is desirable for routine 
work. 


For the ionosphere at UHF and higher frequencies it is a fair approximation 
(with error to be discussed later) to ignore magneto ionic and collision effects 
and take die refractive index in the form 




where 


Consequently 

a 
S 


f^ = local plasma frequency 


2 

= e 

4jr^mcQ 

jcm 

= 80. 614 N 


N, 


_d 

df '~“p' 


n_ ft) ^(fn„) 


)« 


N. 


(3.3-3) 


(3.3-4) 


(3.3-5) 
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But since n = 1+N , n = 1+N "With N jN small in the region of 

5 s H r Sr 

interest. 


N_ 




1+N. 


PS 




-Np 

o 

w 


(3.3-6) 

(3.3-7) 


The difference between these alternative approximations is in any case signifi- 
cantly less than the error in either of them due to ne^ect of the longitudinal 
magneto -ionic term. 


The deviations from these simple approximate relationships may be significant 

in at least two respects: 

2 

1) The 1/f term is subject to higher order corrections the most significant 

-3 

of which adds an f dependence proportional to the main term times the 
ratio of longitudinal gyro frequency to carrier frequency, or as much as 
2% correction at 100 MHz (Ref. 6 ). 

-4 

2) Bending, ignored in the above analysis, introduces an f dependence 
which may amount to a) 0. 8% of the main term at 100 MHz (Ref. 6). 


It is estimated that the simple relation 

0 = -e„ . (3.3-9) 

g P ' ' 

is accurate to wilhin a fraction for frequencies above 100 Mlfe which 

is useful for many purposes. 

It is worthy of note that while ray tracing programs such as REEK are capable of 

-4 

taking into account the f bending term quite accurately, the much more 
significant f term (fa2% at 100 MHz) cannot be treated by any of the known 
operationally practical ray trace programs because that term Is inherently 
anisotropic and birefringent, i. e., allows two, generally coupled, magneto- 
ionic modes, and requires a considerably more complex treatment dependent 
on antenna polarizations among other things . 
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REEK Group Option 

The REEK ray trace program includes an option for group range error computed 
as follows: Recall from section 3,1 that REEK computes the total retardation on 
the Basis 

''p = v + x 

P P . 


where 


'R_ 


= phase retardation error 


e« = phase bending error 
P 

= J* 1 ds - R 
P 

p = the phase path. 


For the group option REEK takes 




where = J(n “-l)ds 

g p 

“ Jlds-R 

g P P 

N_ 


“g-^ ■ = -T+k 


whereas 

g P 

That is the group retardation is computed by integration along the phase path, 
in effect REEK computes 


Pb = 


- 3.14 - 


(3.3-10) 


(3.3-11) 


(3.3-12) 


(3.3-13) 


(3.3-14) 


(3.3-15) 


(3.3-16) 
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Proof of the Group Option Procedure 


This can be justified hueristically on the following basis: 
Starting with the phase range as a function of frequency 


^ lP(f) 

where P(f) the path^ is, by Fermat’s principle, an extremal. 

In other words, considering an arbitrary path deformation, 6{s) (subject to 
the constraint that the deformed patli still passes through transmitter and 
receiver) then it P' = P+X6 where X is an arbitrary multiplier, then the 
integral over the deformed path 


p’ 


must satisfy 



X=0 


J “p* 

P+X6 

0 


for any T satisfying the end constraints (i.e. , 6 = 0 at the ends). 


Now the group range is given exactly from (3.3-17) by the general relation • 
E:jq)anding the path (P(f) about f^ we can write 


m - p(y + (f-y-gSl + o(£-y^ 

whence, denoting f by X, 


0g«o) = / 




ds + 


dX 


■ J 


fyyds 


p(y+x^ +o(x^) 


(3.3-17) 

(3.3-18) 

> (3.3-19) 

(3.3-20) 

(3.3-21) 

(3.3-22) 

0 
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but since the path variation inherently satisfies the end conditions, that 

is, it is an allowed variation of the form 5(s), the path variation of the integral, 
i. e. , the second term above must vanish identically by the extremal principal 
(eq. 3.3-19) and we have just 


Pg(y 


J 


(3.3-23) 


i.e., we have the interesting result that the group range is given by the integral 
of the group refractivity along the phase path. This may be taken as a definition 
of the group path and in this sense it may be said that the group path is identical 
to the phase path in any such arbitrarily inhomogeneous dispersive (but iso- 
tropic) medium. 


Numerical Tests of the Group Option 


As a check on these group -phase relationships a numerical comparison was 
carried out (Ref. 6 ) between the REEK derived group range correction 
as described previously and the theoretically exact procedure of deriving group 
range error by numerical differentiation of REEK derived phase range error as 
a function of frequency. For these tests a hypotlietical Chapman ionosphere 
was assumed leading to an assumed refractivity profile of the form 


N(h) = N ^^^ exp(l -z -exp(-z)) 


where 



16 X 10"^ 

^max 

f2 

^GHz 

z = 

(h~h 

h = 
m 

375 km 

11 

108.333 km 


(3.3-25) 


The REEK phase range error was evaluated at f = .136, .434, and 2.0 GHz 
at various elevation angles. 
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In view of the assumed f refractivity dependence it can be shown by the 

Poincare expansion theorem (Ref. 4, Lefschetz) that any integro “differential 

_2 

function thereof such as e has an expansion in powers of f : 

SpCf) = a^f"^ + a2f“^... 

” •^2 ^4 * * * 


where we define 
A 


^ a f"2 

2 


®4 “ 


In principle then we can solve for the phase ranging error, Sp, in a ray tracing 
program such as REEK at several frequencies, then carry out a numerical fit 
to identify the frequency coefficients a^^ , a .2 or A2 and above, then differ- 
entiate by 3. 3 -2 to give the corresponding derived group error 
e. 


®der 


■Aj - 3B^ 


For precise comparison with REEK group option residts, however, it is 
necessary to take into account that REEK uses the approximation (3.3-15) 

N_ 


while for the f three frequency REEK 


or (3.3-6) for Ng (Ng » " 1+^- 

^ * -2 -2 
results the basic refractivities were scaled in the ratio f , i.e., Np = Kf , 

corresponding to the approximation (3.3-7) (N es»-N ). Expanding (3.3-6) 

®b . P 


N 


^ « -N (1“N„...) 

P P 




•N 

-P P 


So the range error difference on this account is equivalent to a term 
integrated over the path. That is 

"'■g, 




■a Sb 

where for the Chapman profile as defined above this is approximately 
(plane earth approximation) 

2 
P^ 


H. 


pa 




w n: 


max 


HgCscB 


-4 2 

This is an f term as it depends on . 

max 


(3.3-26) 


(3.3-27) 


(3.3-28) 
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The corresponding derived group error is then 



- Aj - 3B^ - 


= -ep-2B^-2H^ 

Which may be compared directly with the group option computed result. The 
comparison is shown in Table 3.3-i in which the columns have the following 
significance 

Sp = REEK phase ranging error, meters 
@ 2B^ - 2 X quadratic fit component of ^ 

(C) 2H^ = 2 X refractivity approximation correction, per eqtn. (3,3-28) 

e = derived group error, eqtn. (3.3-29) 

®der 

REEK group option error, meters. Compare column . 

@ REEK phase retardation and bending terms. In principle these 
should correspond to A2+H^ and respectively. 

The check between tlie multi -frequency derived group error (q) and the REEK 
group option error is veiry close (better than 0. 06%) which would appear to 
. be more than adequate for most applications. The difference though small 
appears to be systematic and numerically significant. Simple modifications 
of the REEK group option have been devised which reduce this error by about 
another order of magnitude in this case but tiae rationale and generality of 
such modifications is not dear and it is questionable whether they should be 
proposed for general use. 


(3.3-29) 


t 
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© 


© 

® 

d) 

® 

© 

Eiev. 

Deg. 

®P 

REEK Phase 

2B4 

= Bending 
Correction 

2 H^ 

= Refractivity 
Correction 

Derived 

Group 

= -@-0-© 

REEK 

Group 

REEK 

Phase 

Retard 

®R 

REEK 

Phase 

Bend 

0.1 

-747.178 

-3.250 

-0.439 

750.867 

751.299 

-749.014 

1.836 

1.0 

-746.213 

-3.143 

-0.439 

749.794 

750.222 

-747.993 

1.780 

2.5 

-741.331 

-2.921 

-0.436 

744.688 

745.098 

-742.992 

1.661 

6.0 

-715.839 

-2.298 

-0.421 

718.558 

718.891 

-717.150 

1.311 

10.0 

-669.932 

-1.626 

-0.394 

671.952 

672.181 

-670.855 

0.923 

15.0 

-603.126 

-1.024 

-0.355 

604.505 • 

604.631 

-603.698 

0.572 

20.0 

-538.671 

-0.656 

-0.317 

539.644 

539.705 

-539.027 

0.356 

30.0 

-435.220 

-0.287 

-0.256 

435.763 

435.774 

-435.368 

0.148 

40.0 

-364. 859 

i -0.134 

-0.214 

365.208 

365.208 

-364.926 

0.066 

50.0 

-318.131 

-0.064 

-0.187 

318.382 

318.380 

-318.161: 

0.031 

60.0 

-287.435 

-0.030 

■ -0.169 

287.634 

287.632 

-287.448^ 

0.014 

70.0 

-268.078 

-0.013 

-0,158 

268.248 

268.246 

-268.083 

0.005 

80.0 

-257,372 

-0.006 

-0.151 

257.529 

257,526 

-257.373 

0.001 

85.0 

-254.794 

-0.003 

-0.150 

254.947 

254.945 

-254.794 

0.000 

87.0 

-254.249 

-0.002 

-0,149 

254.401 

254.390 

-254.249 

0.000 

89.0 

-253.977 

0.000 

-0.149 

254.127 

254.127 

-253.977 

0.000 


In this table, errors are in maters and the following relations apply 


f = 136 MHz = 0.865x10' 

max 

h^^^ = 375 km H = 108.333 km 

max s 


TABLE 3.3-1 

GROUP-PHASE COMPARISONS 
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3.4 Superposition, 


For analytic corrections the tropospheric and ionospheric corrections are 
generally treated as independent. In actuality the incidence an^e for one is 
affected by the refraction error due to the other so that superposition should 
not be expected to hold, i.e, , 


AR 


total 


'^^trop ^^ion ^interaction 


A short study was carried out to investigate the magnitude of the interaction 

error, ®interaction“ done for two different values of ionospheric 

refractivity, -.001 and -.0001. Even with the rather large maximum refrac- 
tivity of -.001, corresponding to say, lOO to 400 MHz frequency, the interaction 

-3 

error is at most of the order of 2 x 10 of the total error at low elevation 
angles as can be seen in Table 3.4-1, where 


®int 


^^tot ^'^^ion ^^trop^ 


The same thing done for range rate errors yields somewhat higher relative 
interaction terms (w3%), apparently resulting from the fact that the tropo- 
spheric and ionospheric errors are of more nearly comparable magnitude in 
range rate. 



(3.4-1) 
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TABLE 3.4-1 


SUPERPOSITION TESTS 


NS = .000313 
RANGE ERRORS 


T-IT AR. 

EL ion 

(DEG) (Meters) 


AR^^„„ AR 
trop ^^tot 

(Meters) (Meters) 


NI(MAX) = -.001 


®int 


RANGE RATE ERRORS 


ARR. 


ARR. 


ARR. 


DIFR ion “ — *trop "*'*‘tot 

(M) (CM/SEC) (CM/SE'C) (CM/SEC) 


®int 

DIFRR 

(CM/SEC) 


0.10 

829.8126 

74.3268 

905. 8979 

1.7584 

00.0000 

00.0000 

00.0000 

0.0000 

0.15 

829.7876 

72. 9693 

904.4672 

1.7102 

0.2871 

15.5816 

16.4216 

0.5529 

0.25 

829.7206 

70.3416 

901.6891 

1.6268 

0. 3852 

15.1146 

15.9795 ' 

0.4797 

0.40 

829.5779 

66.6763 

897.7588 

1.5046 

0. 5494 

14.1084 

15.1284 

0.4706 

0.60 

829.3088 

62.2213 

892.8895 

1 . 3594 

0.7809 

12.9296 

14.1320 

0.4214 

1.00 

828.5021 

54.6030 

884.2277 

1.1226 

1.1815 

11.1572 

12.6854 

0.3467 

1.50 

826.9954 

47.0038 

874.8901 

0.8910 

1.7900 

9.0281 

11.0932 

0.2751 

2.50 

822.3672 

36.2022 

859.1487 

0.5793. 

2. 8147 

6.5692 

9.5735 

0.1895 

4.00 

811.6519 

26.3918 

838.3770 

0.3333 

4.5231 

4.1411 

8.7680 

0.1038 

6.00 

791.3192 

19.0849 

810.5842 

0.1801 

6,8235 

2.4521 

9.3270 

0.0514 

10.00 

736.7049 

12.1096 

748. 8829 

0.0684 

10.1691 

1.298. 

11.4887 

0.0208 

15.00 

658.9756 

8.2814 

667.2837 

0.0266 

13.6325 

0.6714 

14.3112 

0.0073 

20.00 

585.4201 

6.3113 

591.7438 

0.0123 

15. 5324 

0.4160 

15.9515 

0.0030 

30.00 

469.6238 

4.3405 

473.9679 

0.0036 

15.8325 

0.2371 

16.0704 

0.0008 

25.00 

522.0697 

5.1257 

527.2018 

0.0064 

16.0610 

0.3006 

16.3631 

0.0015 

40.00 

392.1420 

3.3828 

395.5261 

0.0013 

14.7871 

0.1828 

14.9703 

0.0004 

50.00 

341.1470 

2.8410 

343.9886 

0.0006 

12.7464 

0.1354 

12.8820 

0.0002 

60.00 

307. 8225 

2. 5143 

310.3370 

0.0002 

10.2844 

0.1008 

10.3853 

0.0001 

80.00 

275.3038 

2.2119 

277.5157 

0.0000 

6.1480 

0.0572 

6.2052 

0.0000 

89.00 

271.6392 

2.1787 

273.8179 

0,0000 

1.7021 

0.0154 

1.7175 

0.0000 



NS = .000313 

NI(MAX) 

= -.0001 




0.10 

82.2783 

74.3268 

156.7865 

0.1814 

0.0000 

00.0000 

00.0000 

0.0000 

0.15 

82.2768 

72.9693 

155.4190 

0.1729 

0.0169 

15.5816 

15,6963 

0.0978 

0.25 

82.2722 

70, 3416 

152.7801 

0.1663 

0.0264 

15.1146 

15.1788 

0.0378 

0.40 

82.2612 

66,6763 

149.0884 

0.1809 

0.0423 

14.1084 

14.2098 

0.0591 

0.60 

82.2389 

62.2213 

144.6006 

0.1404 

0.0647 

12.9296 

13.0250 

0.0306 

1.00 

82.1683 

54.6030 

136.8867 

0.1155 

0.1035 

11.1572 

11.2971 

0.0364 

1.50 

82.0314 

47.0038 

129.1265 

0.0913 

0.1626 

9.0281 

9.2194 

0.0288 

2.50 

81.6000 

36.2022 

117.8598 

0.0576 

0.2624 

6.5692 

6.8521 

0.0205 

4.00 

80.5822 

26.3918 

107.0077 

0.0337 

0.4296 

4.1411 

4.5808 

0.0101 

6.00 

78.6257 

19.0849 

97.7287 

0.0181 

0.6566 

2.4521 

3.1140 

0.0052 

10.00 

73.3061 

12.1096 

85.4225 

0.0069 

0. 9905 

1.2988 

2.2914 

0.0021 

15.00 

65,6604 

8.2814 

73.9445 

0.0027 

1.3409 

0.6714 

2.0131 

0,0007 

20.00 

58,3836 

6.3113 

64.6962 

0.0012 

1.5366 

0.4160 

1.9529 

0.0003 

25.00 

52.0969 

5.1257 

57.2232 

0.0006 

1.5939 

0.3006 

1.8946 

0.0002 

30.00 

46.8824 

4.3405 

51.2233 

0.0004 

1.5741 

0.2371 

1.8113 

0.0001 

40.00 

39.1670 

3.3828 

42.5499 

0.0001 

1.4725 

0.1828 

1.6553 

0.0000 

50.00 

34.0830 

2. 8410 

36.9241 

0.0001 

1.2708 

0.1354 

1.4062 

0.0000 

60,00 

30.7586 

2.5143 

33.2729 

0.0000 

1.0260 

0.1008 

1.1268 

0.0000 

80.00 

27.5131 

2.2119 

29.7250 

0.0000 

0.6136 

0.0572 

0.6708 

0.0000 

89.00 

27.1473 

2.1787 

29.3260 

0.0000 

0.1699 

0.0154 

0.1853 

0.0000 
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4. ANALYnC CORRECTIONS 


This section summarizes the results of a series of studies of the relative accuracy 
and limitations of several proposed and operational analytic and semi-analytic 
refraction corrections. 


4.1 Troposphere 

Table 4.1-1 gives the principal characteristics of the various corrections 
considered here in terms of layer shape assumptions j, measurements utilized^ ap' 
pro^dmations to the integral, and correction quantities computed (R, E, R). 

For the troposphere roughly 90 to 95% of the day-to-day variability is accounted 
for by the variation of Ns, fhe surface value of refractive index. Since this 
is readily measured in terms of pressure, temperature, and humidity it is not 
surprising to find that many of the practical formulations depend on Ns. 

For long range prediction, worldwide statistics of N(z) are available (Ref. 22). 

Reference 18 compared the results of the fo3nuulations listed in Tables 

4.1,2, 4.1.3, 4.1.4 for elevation angle, range, and range rate with the 

corresponding results from REEK ray trace, assuming an exponential model 

"6 

atmosphere with surface refractivity N_ = 313 x 10 and scale height 

s 

H = 6951.25 meters as per the CRPL standard troposphere, 
s 

The results of these comparisons are plotted in Figures 4.1, 4.2, and 4.3. 

Bach of the corrections for either an^e, range or range rate is asymptotically of 
the same form at high elevation an^es, namely, 

AEj .0 = Ng ctn E 

AR^o = HNgCscE 

ARj. = - HN„ E CSC B ctn B 

to s 


TABLE 4.X-1 

TROPOSPHERIC CORRECTIONS 


LAYER INPUT PARAMETERS 

METHOD MODEL (MEASUREMENTS) INTEGRAL R, E, R 


RAYTRACE 







1 

[ 


REEK 

NONE 

N(h) 

RAYTRACE 

R, E 

(Typical of Class) 









SEMI-ANALYTIC 

i 


h 

S 





MOMENTS 

NONE 

M„(h) = r h^Nfli) dh 

ANALYTIC 

R, E. R 




0 


SERIES 



REGRESSION 




dN„ 





NBS STD. 

NONE 

N 

'■^si dh 

TABLEj 

.LOOKUP 

R 


NBS CAPE CANAV. 






1 

R 


WALLOPS 




T, H 



R 


NBS BEAN i GAHOON 



N 

C 



E 


ANALYTIC 






r j 



HOPFIBLD- 

BIQUARTIG 

^^dry 

^ No 
“wet 

ANALYTIC 

R 


DC 

EXPONENTIAL 

N 

s 

ANALYTIC 

R, E, R 

FREEMAN 









NO NAME 



i 






GDAP 









CBAND 







R, E 

NAP-1 







R, E, R 

SAO 







R 


AMS SECOR 






1 

R 


APL TRANET 







R 


NWL TRANET 







• 

R 






( 

1 

1 j 
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TABLE 4.1-2 

RANGE REFRACTION CORRECTION EQUATIONS 


Refraction 

Formulation 

GSFC DC 

GSFC DC 
GSFC Freeman 

GSFC NONAME 
GSFC GDAP 

GSFC NAP-1 

GSFC i SAO lasers 

SAO lasers (after May 1968) 

AMS SECOR 
Wallops C“band 

where AR^.^ = K^IgCscE = 


Range Refraction Correction AR 
AR (meters) = Obs - Corr. 


r 8750 

f .. 1.. \ 

L H 1 

1./+0. 000772 ctn^E / 


AR S750 


H 


AR^ !l - -^ctn"E] 


R. 


r 8432.336 / 1 V 

H ^l + 0.026cscE ) 


f7200/ : 

' 1 "i" 0* 


0045154 csc^E 


® ^ 1 + 0.004csc^E 




°k) 


P 2.238 + 533. 5N_ 

AR T 2—5— 

*- HNg (1 + 10"’^ cos E csc^E) 

^^to["%,^ fl + 0.0236 ctnE ^ 


■] 


A«tc 


r 7600\ 


HAE. secE = first order correction, 
to 
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TABLE 4.1-3 

ELEVATION ANGLE REFRACTION CORRECTION EQUATIONS 


Refraction 

Formulation 


Elevation Angle Refraction Correction AE 
AE (radians) = Obs - Corr, 


GSFG DC 

■ 


GSFC Freeman 



GSFC NONAME 

^to 

f i . A 

I 0.93 + 0.0164ctnE J 

GSFC GDAP 

^toj 

^ A \ 

U' + yi + 0.0045154 csc^E ' 

GSFC NAP-1 

o 

<1 

r 0.000350 f 2 \ 

^ ^s 1 J— J-'J 

® ' 1 + Vi +0.004 csc^E 

Wallops G-band 

^®to 



"I 

J 


where AE..^ = N„ctnE = first order correction, 
to s 
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TABLE 4.1-4 

RANGE RATE REFRACTION CORRECTION EQUATION 


Refraction 

Formulation 

Range Rate Refraction Correction, AR 
AR (meters/sec) = Obs - Corr. 

GSFG DC 

r 8743,25 I 1 

^ ^ ^ (1 + 0.000772 ctn^E)^^^ ^ 

GSFC DC 

/8750'\ 

GSFG Freeman 

A^o " 3?sc^E)] 

S 

GSFC NOHAMB 

r 8432.336 1 1 

^ (1 +0.026 cscE)^ ^ 

GSFC GDAP 

.A {-7200 2 -] 

iO L H y™ -t « J 

VI + 0.00451 54 csc'^E +(1 +0.0045154 esc E) 

GSFC NAP-1 

r 2.7432 2 1 

L HN j ^ 2 J 

Vl + 0.004csc^E +(l + 0.004csc^E) 

APL TRANET 

R „ 

ARj.o -j| sin E [f(E)] 

NWL TRANET 

a6 2.3 
s 

where 


2R 

f(E) = 1 + sinB 

[yH^+2RgH^+R^sin^E’ - R^ sinE 


R (1 + sinE) 

+ (Rg+Hj;) - 

(R^+H^+yH?+2R^H j.+ RgSin^E 


AR^ = -KNgE cscE ctoE = -ARj.^EctnE = ~HAEj,^EcscE = first order 

correction 
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ELEVATION ANGLE (DEGREES) 


H — r- 
SO 


100 


FIGURE 4,1 

RANGE ERROR DUE TO TROPOSPHERIC RBFR^ TION 
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so it is as expected that ail of the various corrections approach a quite accurate 
match to the REEK reference at high d.evation angles , or at least could be 
brou^t into very close agreement ty choice of scale factor^, H. Significant 
differences appear below about 5*^, however. GDAP and NAP-1 (which are 
essentially identical) appear to offer consistently good approximations down to 
the lowest angles. 


References 7, 8, and 19 as sn-nmarized in section 5 of this report derive series 
expansions of the elevation^ angle, and range errors in terms of the successive 
moments of the distribution of refractivity vs height. It turns out that for range 
errors all the way to the horizon and for elevation and range rates errors down to 
about 2° elevation, the moment expansion utilizing the first three (0^, I®*", 2^^) 
moments provides quite a good analytic correction, midway in accuracy between 
those corrections depending only on N_ and the REEK correction based on the 

Q 

entire profile. This is plotted for comparison in chapter 5, Figures 5,7, 5.8, 5.9. 


Hopfield (Ref, 31 ) has proposed a bi-quartic model for the tropospheric refrac- 
tivity; 


N(h) 




where 




surface value of the "dry" terms of refractivity 


77. 6 P 

0 



associated dry teamx sc^e height 



surface value of "wet" term of refractivity 
_ 77.6 (4810 e) 


h. = associated wet terra scale height 

Sff 

Both limited to zero above where they go to zero at hj. and h^. 

a w 

This bi-quartic form fits experimental data quite well over the most important 
part of the troposphere (obviously not around h or h. ) and the resulting 
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straight-line refraction integral can be carried out analytically (Ref. 31 ) 
yielding: 


Ap. = X 

i=l ^ 




t. 

X 


2\^oh ah.%^15 h. 

l I IX i 

+l4(4+l%.4)-4%<4,-h?.) 


2 ^t. ^2 % ^2 ^t.^ 


^i 


Where r^ = radius of earth 

h^ = height of top of layer 


t 

^i 

= 


Xo " 


"3. 

1 


= radius of top of layer 


r siuE 
0 

r cosE 
0 

rt n 1/2 

K - 4 > 

1 


E = elevation an^e of ray at h = 0^ r = r 


■0 


dry term 
wet term 


Actually this closed form is difficult to compute since the individual terms in 
braclcets aire of order j typically 10 , larger than their sum ! While 

this is well within capabEities of some computers in double precision, it would 
be highly desirable to find simplifications of this formula which would make it 
more readEy avaEable on more modest computing equipment. Reference ll 
treats the problem of alternate ways of computing this including two 
very efficient series expansions derived by Yionoulis (Ref. 32 ) and straight- 
forward numerical integration. Best results were obtained with simple 
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numerical integration; using a Simpson formula gave 5-place accuracy with 
20 points for E>2°j 80 points down to 0.6° and 160 points below 0.6°. 

i 

I 

Finally a series cf comparisons was carried out with 85 actual radiosonde 
profiles measured during WICE. 5 methods were compared 


• REEK ra 3 Ttracing (including bending) of the actual profiles (taken as the 
standard of comparison 

• Hopfield*s model 

• The NBS "Standard Sample" regression model 

• The NBS "Gape Canaveral Sample" regression mod<d 


• A special regression model based on T^j derived for the WICE . 
sample. 

A few words of explanation of the latter three models will be in order. The NBS 
regression models are those derived by Thayer and Bean (Ref. 33 ) and are of 

the form 

AR 
where 


Two sets of coefficients were found, one the so-called "standard sample", from 
I a worldwide sample of 77 profiles; the second from 84 profEes all taken at Gape 

i 

Canaveral. 


= A(B) + B(B)AN 

A(E) and B(E) coefficients are functions of elevation angle, 

E, established by least-squares regression fit to a 

sample of actual data and published in Ref. 33 

AN = N^ - tT 
o o 

Nq = ground level measured refractivity 
= average value of N^ over the sample 
AR = predicted range refraction error. 




These regression models were extended for the 85 profEe WICE sample by considering 
various possible regressors of the form 

AR = ^ Ai.{P,T,H,B) 
i 

where H = Relative Humidity 

B = partial pressure of water vapor 
f £(• ) = various tentative fimctional combinations and transformations 
of the basic P, T, H, E measurements. 
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The various combinatious tried and the resulting standard error of the regression 
fit for each case are given in Table 4-5. 

TABLE 4-5 


STANDARD ERROR OP REGRESSION FIT, METERS 


Independent 
. Variables 

^i> 

Dependent Variables 

AR(E = 12.173°) 

AR(E= 90°) 

AR(E= 90°) 

.032 

— 

T 

.216 

.046 

E 

.162 

.035 

N(P, T, H) 

.149 

.032 

E, P 

.154 

.035 

B, P, H 

.147 

.031 

T, H 

.159 

.035 

T, H, P 

.148 

.032 

P/T, E/T^ 

.148 

.032 


The particular case where AR(E = 90°) was used as the regressor (for lower 
elevation angle) is of some interest in indicating the possible residual error at 
other than vertical incidence if the integrated refractivity content (equivalent to 
AR(E = 90°)) could be measured by some other means; die indication is clearly 
very encouraging to such a development. The upshot of these studies was to 
indicate that given the several ground level measurements (P, Tj H or E) about 
the best and simplest predictor is N(P,T,H) alone^ i.e. , just the form of the 
NBS model, and this was the form used in the comparison. The coefficients 
resulting from the three regressions are compared in Table 4-6. 
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TABLE 4-6 


REGRESSION COEFFICIENTS AT B = 12.173° 




®(m/N) 

NBS Standard Sample 

3.8929 

.02149 

NBS Cape Canaveral Sample 

6.6595 

.01377 

WICE Sarajie 

7. 3483 

.01100 


AR=^ = A + BAN 

Finally the results of these various regression models along with the Hopfield 
bi-quartic model were compared with REEK ray tracings for the 85 WICE 
radiosonde samples. For this purpose the scale heights of the bi-quartic model 
were taken as hj„, = 41.17 kra and h ^ = 12.0 km as recommended by 
Hopfield for the latitude of Wallops, The wet, dry^ and total refractivity terms 
in all cases were taken from the ground level readings of the radiosonde. 

The results are summarized in Table 4-7 for die case of E = 12°. 


TABLE 4-7 


COMPARISON OF TROPOSPHERIC REFRACTION CORRECTORS 
VS REEK RAY TRACINGS ON 85 WICE PROFILES 


Predictor 


Std.Dev. of Prediction Errors (Meters) 
E = 120 


Hopfield 

MBS "Standard Sample" Coefficients 
NBS "Cape Canaveral Sample" Coefficients 
Special Wallops Regression on N(P,T,H) 


.223 

.281 

.258 

,149 


It should be pointed out that the Special Wallops Regression is undoubtedly 
optimistic in its estimate of its own prediction error since, among other things, 
it was developed and tested on only a three -month segment of data (April, May, 
June 1968), In comparison with the NBS predictions it will be seen that the 
Hopfield formula is generally comparable, but slightly better in accuracy. 
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4.2 Analytic Correctioas for the Ionosphere 

Reference 20 presented a comparison of several commonly used anal 3 rtic 
corrections for the ionosphere with REEK ray trace taken as a standard. 

The external characteristics of these algorithms are given in Table 4-8 in 
terms of layer model assumptioas ^ input parameters, integral approximations 
and outputs. 

Tables 4.9 - 4.11 give the actual formulations of the analytic corrections and the 
moment correction is discussed in section 5. For the comparisons shown in 
Figures 4.4 - 4.6 a Chapman -type ionosphere was assumed with = 10.67 xlO 

h^ = 364 km and = 104.667 km, typical of daytime near solar maximum at 
f = 2GHz. The GEOVAP and Freeman formulations are not useful at low elevation 
angles and are in fact relatively poor at angles as high as 40° - S0°. All the 
other formulations listed in Tables 4,9- 4.11 except NAP-3 act simRarly at 
low angles in R, E and R and what differences do appear can largely be ascribed 
to the necessarily somewhat arbitrary choices involved in fitting the parameters 
of the various models to the particular Chapman model used as the basis of 
comparison. 

In other words, a slightly different choice of h^ and parameters would 
in these cases yield vary close agreement with the REEK ray trace results. 

The ionosphere is significantly different than the troposphere in that the minimum 
elevation angle in the ionosphere is of the order of 18° - 20°even for horizontal 
takeoff from the earth (E = 0), consequently the very low angle problem is never 
encountered in the ionosphere for the cases with which we are concerned. 

Consequently the main difference between approaches is related simply to how 
much information about the actual ionosphere is actually input to the algorithm, 
or how well the modelled ionosphere matches the real. 


Moment series comparisons with idealized and actual ionospheres are given in 
section 5. The match is very close as would be expected from the fact that the 
moments themselves approach a complete description of the actual layer and 
the very low angle convergence question is moot for the ionosphere for the reasons 
discussed above. 
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TABLE 4. 8 

IONOSPHERIC CORRECTION ALGORITHMS 


METHOD 


LAYER 

MODEL 


INPUT 

PARAMETERS 

(MEASUREMENTS) 


INTEGRALS R, E, R 


RAYTPJICB 

REEK 

SEMI-ANALYTIC 


MOMENTS 

ANALYTIC 
DC 
DODS 
Freeman 
GEOVAP 
GPRO (GDAP) 
NAP-3 


NONE 


NONE 


CHAPMAN 


PARABOLIC 
Parabolic BottomI 
Exponential Top 


N<h) 


MOMENTS 


^maXj ^max 


y_= bottom scale ht 
*'m 

k= topside scale 

height of max 
m ^ 

^max 


FULL RAYTRACE 


ANALYTIC 

SERIES 


ANALYTIC 


R, E 


R, Ej R 


R, E, R 
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TABLE 4.9 

RANGE REFRACTION CORRECTION FORMUI.AS 


Refraction Range Refraction Correction AR 

Formulation AR (meters) = R ~ R^ 



- -Ni^He^cscE = -M^cscE 
10 im o 


4,16 - 
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TABLE 4.10 


RANGE RATE REFRACTION CORRECTION FORMULAS 


Refraction 

Formulation 


Range Rate Refraction Correction AR 
AR (meters/sec) = R ~ R„ 


GSFC DC i DODS 

GSFC DC i DODS 
(E>10®) 

GSFC Freeman 


GSFC GPRO (GDAP) 


-3/2 


4 . / ^e \ . 3 

^ ^^io( Rg+h j ® ^ 

\ e m / \ e tn 


4 - 

5 10 


r / o X/H+h^ 

\ sm E /V e 


1 - 4 P^ 

■8N. H„E, ctnE, cscB, 1 -j s s 5-7? 

^ Z ^ i (l+4s2ctn2Ep^/^ 

[1 +(1 + 4^2 ctn^ 


NAP-3 


/ K V 

E (AR) ' I sin B cos E 


R„ cos E 
e 


GEOVAP 


^^io 


AR. = N. He EctnEcscE = M_ • 
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TABLE 4. a 

ELEVATION REFRACTION CORRECTION FORMULAS 

Refractioa Blevatkm Angle Refraction Correction AE 

Formulation AE (radians) = E — E 





ERROR (METERS) 





(MlULiaEGREES) 




RANGE RATE REFRACTION ERROR (CM/S EO) 
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DC/3X)DS(>10°) 

f ^ ®io 


AVERAGE 

364,000 km 

-Njjj^=10.G7xl0“® 

H “ 104.667 km 
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GEOVAP/AR 


m (GPRO) 


NAP-3. 


m (GPRO). 


m (GPRO) 


GPRO''^ 
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FREEMAN - 


iREEK (BENT) 


^~REEK (CHAPMAN) 
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5. MOMENT EXPANSIONS 
5,1 Introduction 

The refractive range errors at a given frequency - (to a zero order approx- 
imation valid for plane earth) - are directly proportional to the total columnar 
r^ractivity content or zeroth order moment of the density profile and a simple 
function^, esc E^ of the devation angle. For spherical earth, a first-order 
correction is available which depends only on the eiSective height, or first- 
order moment of the electron density profile. 

These considerations lead one to consider the possibility of an expansion of tb,e 
range error correction in terms of successive moments of the profile dis- 
tribution, and in which the coefficients are elevation an^e functions. If 
such a formulation is capable of providing sufficient accuracy in a reasonably 
small number of terms, it would have significant computational advantages 
over straightforward raytracing for the case where corrections at a large 
number of points (i. e. , elevation angles) must be computed for the same 
ionosphere. 

Freeman (Ref. 26) made a start in this direction for ionospheric range errors, 
but found relatively slow convergence. His expansion was in terms of 
momencs about the ground level, h = 0. It can be reasoned that the expansion 
would necessarily be more rapidly convergent if moments were taken about 
a point in the ionosphere, such as the centroid. For example, it is clear, 
to the extent that the ionospheric density profile can be approximated by a 
thin shell at its centroid, that a single term in the moment expansion (the 
zeroth moment) would give an ejsact answer for the case of moments about 
the centroid, since all higher order moments are already zero for moments 
about the centroid, but not for moments about any other point. 

This section summarizes further residts on the moment expansion first 
derived under this project and reported in Refs, 7,8, and I9, 
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5.2 Moment Expansion Derivation 

Based on the straight path assumption^, discussed in section 3.2, the range 
error is approximated simjiy as 

AR = J* (n-l)ds 
strai^t 
path 

= J’aec$^N(h)dh . 


On a spherical earth (See Figure 5,i) this becomes 

*. ^ • 

AR = J sec(J(h))N(h)dh 
o 

where N(h) = n(h) - 1 

$ (h) is the angle the straight line ray makes with the local 
vertical as it passes through height h. 

From the spherical geometry shown in Figure 5il , sec§ (h) is given by 

> 2 -, - 1/2 


secj (k) = [l - ] 


Now, following the motivation discussed previously, we expand sec ($ (h)) 
in a .Taylor series in h about some arbitrary reference height, b^, which 
will normally be chosen somewhere near the center of the refractive layer, 
thus 

sec(«0i)) = sec(|(h^)) + (h-h^) [~sec(§(h))]... 


^ (h-h^) 


m 


^m 


m! 


m=0 


dh^ 


m 


(see® (h)) 


where the derivatives are to be evaluated at h = h , 

c 


(5.2-1) 


(5.2-2) 


(5.2-3) 


(5.2-4) 


Note that the geometrical factor is independent of h. This means diat when 
we carry out tlie integral, equation (5.2-2), interchanging integral and 



Communications Research Laboratories 


summation, the geometrical factors come out of the integrals and we have 



where 


m 


am 


ml 


sec(§ (h)) 


are geometrical factors 


h=h 


depending only on cp(h = h ), i.e, , on cp and h . 

^ . c o c 

**8 

= ; (h-h^)”NCi)dh 

0 

= mth moment about of refractivity distribution. 

The moments adequately characterize the ionosphere or 
troposphere for refraction correction purposes, as will 
be shown. 


5.3 Moment Expansion for Angle 


The corresponding moment series for elevation angle errors can be derived 
from that for range from the general relation between range and angle errors 
derived in Ref. 13, namely 


AE = N^ctn E + (i + esc E ~ AR (E, h) 


where 


a = radius of earth 
R = range to target 


~J({ 


2 2 

Ja+hj.) - (a cos E) - a sin B 
hj. = target height 
Njj = ground level refractivity. 

Using eq. (5.2-5) in eq. (5.3-1) gives the moment series for angle 


AE w N ctn E + 




cscB y G;^(B)M^C^) 
m=0 


where g;„<B) = 

The first few G’ and G are given in Table 5.1. 


(5.2-5) 

5.2-6) 

(5.2-7) 


(5.3-1) 

(5.3-2) 
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TABLE 5.1 


GEOMETRICAL FACTORS 

9 

C = (acosE) 
r = a + 


HI 


m 




1 --^ 
r 2 


rl /2 


-CG^/r^ 


3/2CG^/r* 


^GC-SG^+G^) 

-^(7G^ -3G^) 
^ o o' 


3C 


^(- 2 iG^^+l 4 G^-G^) 
7 o o o 


8 r 


-^(33G^^ -30G^^ + 5G^) 
o o o' 


G* 

m 


, - CG^ taa S/r^ 

G* p 3CG^ 

lfor2Q2+ 1 


GI 


J-r [2<-5 G*+g2)+ 4(-35G«+5G^) 


G* 


|-r[2(7G« -3G^)+ 4(630® -21G^); 


^-|.[2(-21G® + 14g“-G^) 


+ (“231 + 126 G^ “ 7G® 

O 0 0 

r 


G’ 


5^4[2(330f -30G® + 5G«) 


+ (429gJ;^ - 330gJ;° + 45G? 

^*0 0 0 
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5.4 Moment Expattsion for Range Rate 

The momeat series for range rate is given simply by differentiating 
eq, (5.2-5) 

AR = E ^ G;^(E;)M^(h^) (5 

m=0 

This assumes that the moments themselves are not changing, i. e., that 
the target is either well above the ionosphere or traveling at a constant 
level in the ionosphere and that the horizontal variations are negligible. 

This assunption also underlies die other refraction correction analytic 
developments discussed in this report. Note that the series here is iden- 
tical to that for the elevation angle so the convergence behaviour is identical. 


5.5 Convergence 


Convergence of the moment series was considered in Ref. 8 for 
the specific case of tropospheric elevation angle where convergence dif- 
ficulties had been noted. The general result gives a sufficient condition 
for convergence as 


a+h 


cosE 

For example: 
Ionosphere 
assuming 


a 




a = 
h = 


6380 km 
375 km 


Hg = 108 km 
we find cosE < 1.05 

which is of course always satisfied, i.e. , the series should always 
converge for the ionosphere. 
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Troposphere 
using a = 


•we find E > 


6380 km 
0 

6.95 km 

l.B9° for convergence. 


We shall show examples indicating diat the seines doeSj in feet, diverge 
below about this value. 


5.6 Tests of the Moment Expansion for Range 
5,6.1 Ionosphere 

Three types of test of the moment expansion have been carried out. Each 
involves comparison of the moment expansion against raytrace calculations 
by the REEK program for various ionospheric density profiles at the SECOR 
effective frequency of 434.2696 MHz. 


Figure 5.3 shows the residual range error as a function of the order for 
expansion at the centroid compared to REEK ray timees through the first 4 of 
the 17 WICE ionospheric profiles described in section 6. This makes it clear 
that the convergence is far from uniform and that the expansion should be stopped 
on even order with N = 2 adequate. 


in the first test (Ref. 7 ) a model ionosphere was postulated as a modified 
Chapman layer of reasonably representative thickness, height, and density. 

This was used mainly to explore the convergence properties of the range 
error expansion as functions of the expansion height, h^, and order N. Figure 

5,2 shows the residual error of the partial sum ) G_ (E )M_ to terms of 

^ m' o' m 

order N. for various N. as a function of the expansion center h . Clearly 

t L C 

the convergence is best for h at or slightly above the centroid height. 


- 5.7 - 



Couvergence Error, inetors 


Communications Research Laboratories 


ipppj-wwpiniiinpppp — Hp»p»p— PfpM—p wiip—— ■ ■ ■ ■■■Up w a p pppppppp 

<wwpiiip» Pit il#5p»p*p< iWwi*g !M>—fpa>M»»— B— jPjP 

|u»pnr»H|U^u»:S»i^^sgH^ 


■■ «M>PPB nMHi ■ MHt K p P ■ 




:»ra 


■hpbppbi 

IppaPPPw 

!■■■««■■ 


■pgpppp 


. ;sasff 


:3s»sa: 



:aama 







CHAPMAN IONOSPHERE 

N 

max 

= 

200 • 10"^ 

^max 

s 

375 km 

Hs 

=5 

108 km 

H 

sat 

:= 

1333 km 

H 

centroid 


437 km 



Hi B«i9*P ■*■*!■ ■■•«• 


lasHiiMPM* 




S\ 


:3Ks 


SmSppi 

mM 


aaaaffi! 




I ■■■■« BPrtP *«< 


pauiuup 





iPPWip pCT iiPi 
■■nwiPMMl 
■ ■VHPPPiWBB 
■sBBBfMtlPSPB 


■ ■BBBBBIMaPi 

■ BPBMiBPBi 

BPNIIIMIMB 

PBBMPJgg 




■BBHBiSaiBi 

■BBBBBBglggP 


, Jbs 

JBBBa 

■aPWiiNai* 

Ibpppbbbbmp 

iBPBBBBBBiai 

■ BBBBBBHU|r 
■ttitiiPBpniir; 
■lIPWItBSMUk 

■ HBPPBBPMiE 





FIGURE 5.2 

CONVERGENCE ERROR OF THE IONOSPHERIC 
MOMENT EXPANSION FOR RANGE VS EXPAN- 
SION CENTER h FOR VARIOUS ORDERS OF 
EXPANSION, N^ 

ELEVATION ANGLE, E = 10° 




Expansion Center, h km 

C ’ 


ORIGMAL PAGE 










Communicatious Research Laboratories 


The remaiaiQg test cases detailed in Ref. 7 were based on real profile data for 
16 of the 17 WICB ionospheric data sets. These were analyzed in the same way 
with respect to partial sum convergence. The results were completely in con- 
sonance with the above findings for tlie Chapman layer and need not be repeated 
here. 

Figures 5.4_, 5.5 and 5.6 show the moment series range, elevation, and 
arange-rate corrections for the same standard ionosphere used in Ref. 20 
for comparison witli otlier commonly used analytic approximations and the 
REEK raytrace as a standard. It can be concluded that the moment series 
to N = 2 provides a convenient and more accurate analytic approximate 
correction than any of the other commonly used analytic forms for ionospheric 
range, elevation, or range-rate tested here. 


5.6.2 Troposphere 

Tests of the moment expansion for tropospheric refraction error were carried 
out on an exponential tropospheric model using REEK Raytrace as a reference. 

I 

The results shown in Table 5.2 as a function of expansion order confirm the 
conclusion previously derived that the series for the troposphere becomes 
divergent for larger orders below about B = 2.5°. The partial sum to 2 terms, 
hawever, provides a simple and reasonaldy accurate representation. 

As additional tests the series to Nj, = 2 for range, elevation and range -rate 
were calculated for the same standard troposphere used as a basis for comparing 
various analytic approximations in Ref. 18. The comparative results are plotted 
6t Figures 5.7, 5.8, and 5.9. Agreement with the REEK ray tr.. es is better 
fiarthe moment expansion corrections than for the other commonly used analytic 
approximate correction forms tested here, at least down to about 2° elevation 
angle. Some of the analytic expressions are better than ths moment expansions 
below 2° for tropospheric refraction corrections for elevation and range rate. 
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TABLE 5.2 


TROPOSPHERIC RANGE ERROR 
ARt’ (ill) 


ELEV. 


MOMENT EXPANSIONS TO 

ORDER: 


REEK 

Degs. 

0 

1 

2 

3 

4 

5 


0.1 

46.609 

46.609 

64.023 

35.035 

149.033 

351.920 

V 

0.5 

45.842 

45.842 

61.868 

36.062 

134.230 

283.065 

63.6367 

1.0 

43.684 

43.684 

56.273 

37.865 

101.450 

143.975 

54.1222 

2,0 

37.359 

37.359 

43.113 

36.961 

52.501 

8.640 

40.7921 

4.0 

25.946 

25. 946 

26.872 

26.395 

26.976 

26.185 

26.3309 

5.0 

22.028 

22.028 

22.435 

22.284 

22.417 

22.287 

22,1553 

7.0 

16.689 

16.689 

16.789 

16.768 

16.779 

16.773 

— 

10.0 

12.112 

12.112 

12.132 

12.130 

12.131 

12.131 

12.1028 

20.0 

6.310 

6.310 

6.310 

6.310 

6.310 

6.310 

6.3105 

40.0 

3.380 

3.380 

3.380 

3.380 

3.380 

3.380 

3.3811 

80.0 

2.209 

, 2.209 

2.209 

2.209 

2.209 

2.209 

2.2092 

90.0 

2.176 

2.176 

2.176 

2.176 

2.176 

2.176 

2.1756 
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N_ = 313 . i0“‘ 


H„ = H„ = 6951 m 
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6. PROFILE DETERMINATIONS 


Because of the predominant role of ionospheric errors in the systems under 
consideration, a considerable effort was devoted to a study of how wdd the 
ionospheric profile could be predicted or measured by availabie techniques, 
including bottomside and topside sounders and predictions. This section 
describes the techniques used for these predictions and their comparisons. 


6.1 Jacksou*s Composite Profiles 

The principal source of ionospheric prof Re data was that provided by John E. 
Jackson of GSFC for 16 of the 17 Alouette coincidences identified in Table 6.2a). 
These consisted of composite bottomside and topside profiles, based on true 
height reductions of bottomside ionosonde and topside Alouette data. In some 
cases these involved both Alouette 1 and Alouette 2 passes and in some cases 
both Ottowa and Grand Bahama Island bottomside soundings which provided a 
basis for estimation and incorporation of Latitudeinal gradients. Alouette 2 data 
by itself was of limited usefulness because of its low altitude («i530 km), only 
slightly above some adjustment was necessary to the topside 

heights to account for the observations discussed elsewhere (Ref. 24) that 
topside ionosonde reduction apparently tends normally to estimate true heights 
that are too low and bottomside-topside composite profiles meshed straight- 
forwardly tend to predict a total ionospheric content (or effective thickness) 

2% to 15%, averaging about 5% too low. The details of how these various 
complications were handled are recorded in a series of GSFC memoranda 
and documents (Refs. 24 and 25), Suffice it to day here that these "Jackson 
composite" profiles represent what is believed to be the best product, of 
skilled judgement based on available topside and bottomside soundings. 
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6.2 ESSA Bottomside Profiles 

ij 

'i 

For comparison v/ith these "Jaclcson composites" ESSA was asked to do an |i 

independent reduction of the associated bottoms ide-only profiles for six days - 

for which good ionosonde data was available. ESSA ftirtiier provided estimates * 

of the entire ionospheric profile based on extrapolations of the bottom- I, 

side-only data using a modified Chapman layer fit. 

i’ 

I! 

i! 

6.3 Predictions from Radio Propagation Predictions !; 


A method that has been in use at GSFG for some time utilizes radio pro- 
pagation predictions as an indirect basis for estimating ionization density 
profiles . The origin of the method is unknown but it is disci issed and refined 
in terms of its numerical parameters in Ref. 26 ^ which is the form adopted 
for this sutdy. As a starting point the quantities 


and 



M(3000)F2 = 
or 

EJF(0)F2 = 


and 


F 2 layer ordinary critical frequency 

- ' * * . 

3000 km MUF factor 

MUF (0)F2 = Estimated Junction Frequency or "Standard" 
Maximum Usable Frequency reflected 
from F 2 layer at vertical incidence 


EJF (4000)F2 ” MUF (4000) F 2 = same thing at 4000 km range 

¥ 

are obtainable from predictions, e.g. , Ref. 27, either numerically or graphically. 
In the latter case, the EJF quantities are normally given in which case these are 
transformed to the former two quantities by (Ref. 27) 


V 2 

and 

M(3000)F„ 


% 


= MUF(0)F2 - -~ 

EJF (4000) F 2 
iaxfoF2 


where % = local magneto -gyro frequency 


eB 

2jtm 


« 1.5 MHz at Wallops Island. 
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In terms of these available predicted quantities j then, the layer refractivity 
is modelled by a modified Chapman shape 

N(h) = N^exp(l-z-e“^) 

h-h 

■where z = — jr 

^s 

with parameters related to the predicted quantities by 



h^ = 1393. 1 exp (-0.5014 xM (3000) F^) 

(= height of layer maximum) 

H, = |[30+.2(h„-200)] 

h^-50 

3 

For future reference these equations will be called the Freeman prediction. 


Note that under this model the total refractivity integral or zeroeth order 
moment is given by 

JXfdh = 

= N„H e 
in s 

where e = 2.71828 ... 


Two different sets of basic ionospheric predictions were compared under 
Freeman's method, the first using the "Basic Radio Propagation Predictions 
Three Months in Advance", published for the period of the tests by CRPL 
(now ESSA) and the second based on CCIR repojtt #340, "Atlas of Ionospheric 
Characteristics", (Oslo 1966), using the relevant sunspot number R = 112. 

Both pro'vide the same type of data but with slightly different numerical results. 
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6.4 BSSA Predictions 

For some time BSSA has provided to GSFC monthly predictions of maximum 
refractivity and h^ wbrnh are then used in a form equivalent to the 
above equations for minitrack corrections. The predictions are valid for the 
month but for simplicity are given at only 4 times during the day which are 
then interpolated linearly to the time of interest. 


6.5 Comparisons in Terms of Moments 


A convenient mode for comparison was available in terms of the moment expansion 
technique devd.oped in section 5 whereby 


N 

= 1 

m=0 

where M = m^ moment of refractivity profile distribution in height 
. \ 

= J“ (h-h^)™N(Ii)dh 

o 


h_ = satellite height^ taken as infinite for these comparisons 
N = expansion order. 


From section 5 it was found that the expansion to N = 2 suffices so that 
only the N = 0, 2 teams need be considered. 


The moment expansion has the merit for these purposes of explicitly 
separating the angular dependence from the profile shape dependence in the 

G^(Eo) and terms, so makes possible a comparison of fundamental 

layer parameters independent of particular ray geometsry, in terms of the 
first three moments, M^, and Mg. Note that the dominant term, M^, 
corresponds simply to the range error on a vertical ray or total refractivity 
integral. The first three geometrical coSBficients, Gq(E), Gj^(E), Gg(E) 
are plotted in Figure 6.1. 
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The moment series comparison was carried out in detail for the "Jackson 
Composite" and "ESSA Bottomside" profiles in Ref. 28 for 6 days. 

The ESSA profiles were given at 3 to 5 15 minute intervals j about the reference 
time taken as time of point-of -closest- approach. 

The first 3 moments and were found by numerical integration 

for each adjacent time for each day, giving 

n,t,d 

n = order = 0,1,2. 
t = record index for day = 1,2,... T^ 

T^ = 3 to 5 
d = day = 1,..,6 


For the purpose of this analysis the 3 to 5 data within +45 min of TCA for each 
day were averaged giving 



which were then compared with the corresponding moments for the Jackson 
composite profiles, denoted M^, yielding the differences 


DD 


n,d 


= - mJ, . 

n,d n,d 


The mean and standard deviation statistics of DD^ ^ were computed yielding 
the data of Table 6 . 1 . - 


It is notable from the Table that the sample bias term is only of the order 
of 1/3 of the estimated standard deviation, in other words, the discrepancies 
between the two types of data do not appear to be significantly biased. 


The significance of the various orders of moments in terms of the total range 
error discrepancy is shown in Figure 6,2 where the standard deviations above 
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TABLE 6.1 

COMPARISON OF BSSA BOTTOMSIDB PROFILE 
EXTRAPOLATIONS AND JACKSON COMPOSITB 
PROFILES IN TERMS OF MOMENTS 
f == 434 MHz 


ORDER OF MOMENT 






a = 0 

1 

2 



(meters) 

(meters)^ 

(meters)^ 

DDn 

= Mean 

-1.79 

8.62 E4 

1.53 ElO 


= Std.Dev 

3.075 

2.34 E5 

3. 97 EiO 


= RMS 

3.56 

2.49 E5 

4.25 BIO 


t 
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are multiplied with their corresponding coefficients and plotted as the terms 

OqGq(E)^ Gj^(E), ^2 G2(H). 

It is immediately evident from Figure 6.2 that the zero order term predominates 
by more than an order of magnitude and accounts for better than 99% of the total 
variance even in the worst case of low elevation angles and fully correlated 
moment discrepancies. The analysis of Ref. 28 considered the actual correlation 
effects between the various orders of moments by computing the full moment 
discrepancy covariance matrix and propagating the resulting covariance in the 
correct manner through the moment series summation to compute the total 
ranee variance; the results reconfirm the above general conclusion as to tlie 
dominance of the zero order term. 

The first and second order moments add significantly to the accuracy of the 
zero order term for actual data correction. For the purpose of profile accuracy 
comparisons, however, the above analysis indicates that essentially all the 
significant difference is in the zeroeth order term and further comparisons can 
be made most efficiently in terms of alone, independent of elevation angle. 

This comparison is carried out in Table 6.2 for the 17 Alouette-GEOS coinci- 
dence passes, and for the four predictions discussed above. In addition, as a 
point of comparison. Table 6.2 includes the actual ionospheric total content 
measured by the S3COR two -frequency data at the point of GEOS closest 
approach and adjusted from measured slant path to vertical path by the geom- 
etrical coefficient G^(E) evaluated at the elevation angle corresponding to the 
point of closest approach (50 degrees or greater in all cases) and using, in 
this case, expansion center h^ = 350 km. 

The rms differences of the various predictors (columns) in Table 6.2 a) are 
shown in Table 6,2 b). The rms discrepancies lie generally in the range 
4 to 11 meters or roughly 15 to 35% of the total error. It is tempting to utilize 
the matrix of rms discrepancies as a basis of estimation of the absolute errors 
of each of the columns in a modified Grubbs variance analysis. This is probably 
not legitimate because of correlations between the errors of various predictors; 
for example, columns 2 and 4 (Freeman predictions based on different ionospheric 
radio predictions) should be expected to be higbjy correlated. 
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TABLE 6.2a) 


4 

4 

4 

4 

4 

4 

5 
5 
5 

5 

6 
6 
6 
6 
6 
6 
6 


COMPARISON OF IONOSPHERIC PROFILE PREDICTORS IN 
TERMS OF M^ (TOTAL REFRACTIVITY INTEGRAL) 
(METERS) AT F = 434 MHz 



1 

2 

3 

4 

5 


SECOR 

ICOR 

(Meas urement) 



CO 

c 

Q 

bo CVl 

CO 

a 

GMT 

>\ • c 
Q X 2 

F reeman 
Method 
Using 3 Mo, 
Predictions 

Modified 

Jackson 

Composites 

Freeman’s 
Method Usir 
CCIR Predic 
For R = 11 

ESSA 

Bottoms ide 
Profile 
Extra polatic 

03 0143 

49.14 

32.51 

43.30 

35.80 

47.49 

05 0221 

33.49 

29.60 

29.75 

33.00 

28.69 

10 0207 

36.12 

30.21 

33.21 

33.60 

31.61 

12 0245 

27.44 

28.55 

28.86 

32.10 

30.37 

17 0231 

35.27 

29.16 

32.31 

32.60 

38.69 

22 1743 

##.## 

58.13 

##.## 

73.30 

##.## 

24 0320 

51.40 

26.95 

46.58 

29.20 

47.22 

25 0339 

31.47 

26.09 

27.26 

27.97 

##.## 

29 0306 

42.08 

##.## 

33.11 

30.11 

##.## 

30 0325 

36.45 

26.70 

28.64 

28.88 

##.## 

04 0312 

34.06 

26.12 

27.60 

27.90 

##.## 

05 0331 

37.47 

25.04 

35.19 

27.63 

##.## 

11 1715 

29.03 

##.## 

26.78 

33.83 

##.## 

13 1753 

33.34 

##.## 

30.53 

35.91 

##.## 

18 1739 


##.## 

42.15 

35.15 

##.## 

21 0458 

23.07 

##.## 

19.30 

21.55 

##.## 

25 0425 

33.09 

##.## 

28.73 

24.76 

##.## 


TABLE 6.b) 

RMS DIFFERENCE MATRIX (METERS) 


PREDICTION # 



1 

2 

3 

4 

5 

1 

0 

11.42 

4.78 

8.77 

3.74 

2 

11.42 

0 

7.94 

5.32 

11.04 

3 

4.78 

7.94 

0 

6.09 

3.28 

4 

8.77 

5.32 

6.09 

0 

9.34 

5 

3.74 

11.04 

3.28 

9.34 

0 
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Probably a more meaningful interpretation is to talce the SECOR as a reference 
standard on the basis that its error is believed tcTEe^ell under the 3 or 4 meters 
3 nns discrepancy between SECOR and the best of the predictors. Interpreted in 
this way, rough upper limits on the absolute error of the various predictions 


are as follows 

Am 

A 

rms 

Std. 

Dev. 

^A 

RMS 

Profile Measurements 

meters 

meters 

meters 

% of in 

Jackson Composite 

-4.1 

4.8 

2.5 

13 

ESSA Bottomside 
Extrapolation 

(only 6 cases) 

-1.5 

3.7 

3.4 

10 

Profile Predictions 
CCIR Report 340 

-5.2 

8.8 

7.1 

24 

3 Mos . Advance 

Predictions 

-9.1 

11.4 

6.9 

32 


The assessment of the ESSA bottomside profile extrapolations is of questionable 
validity as it only includes 6 cases. The Profile measurements clearly provide 
a better basis than the predictions with errors something less than 10-15%. The 
Freeman predictions based on the CCIR #340 radio predictions show about 24% 
residual error, remarkably consistent with the original claim of Freeman (Ref. 26) 
of 25% but notably on a completely independent data set. It may be marginally 
significant that all the mean biases are of a negative sense. 


6,6 Geographical Gradients 

The visibility period of the GEOS satellite during WICE typically was such that 
the line -of -sight ray swept over some 20 degrees of latitude through the ionos- 
phere. Naturally, the ionospheric changes over this distance are often signi- 
ficant, occasionally as much as 2:1 in terms of total electron content or M^. 

Considerable effort had been expended on defining the Jackson composite topside- 
bottoms ide ionospheric profiles at Wallops Island with considerable success 
when compared to SECOR measurements at the point -of -closest approach to Wallops 
as reported in the previous section. An investigation was then undertaken of various 
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proposed methods of scaling the Wallops predicted profile for its geographic 
variation based on either predicted or observed (bottomside sounder) geo- 
graphical variations. 


Define 

N^(W,h) 


p(e,h) 

NP(p,h) 


(p) 


(p) 

max'^' 


h. 


Jackson composite refractivity at Wallops at Height h, 
applicable at time of closest approach. 

geographical position ^at,long) of height h along ray at elevation e, 

predicted (Freeman -CRPL Radio predictions) refractivity at 
position p, height h. 

predicted maximum refractivity at point p - proportional to 

(V/- 

2 

measured max refractivity at proportional to • 

reference height in the ionosphere^ taken as 425 km. 


Definition of scaling methods; 




r wP(p(e,h),h) 1 

^ N^(Wjh) ^ 

(Scaling by predicted ratio of 
electron densities at height h) 

Ng(h) = 

NVv,h) 


(Scaling by predicted ratio of 


(£ 0 ^ 2 ^^ at p(e,h)) 

Ny(h) = 


(p(e,h))_, 

nP (w) 
max' * 

(Scaling by predicted ratio of 
^^ 0 ^ 2 ^^ at p(e,h^) ) 

Ng(h) = 

N^(W,h)[l 

. P-(p(e,h„)) -L„) 

(L ) 
max' w' 

-* -1 

- 


(Scaling by Latitudinal Interpolation 




2 

between observed ) 

Note that in N„ 

and Ng the scaling factors 

are not functions of height so that 

the ray tracing can be done on N^(W^h) and the result scaled by the factors 


Fj or Fg in brackets. For #8, the scaling was based on linear Latitude inter- 
polation in terms of refractivity (or f^) between the Ottawa and Grand Bahama 


^0^2 
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Each of the corrections; 

R = SECOR measurement 
o 

R = Ray trace on local composite profile ~ not scaled geographically 
c 


Rg = Ray trace on Ng 
®7 = 

Eg = R^xFg 


were computed at the southerly and northerly points of 25 elevation and point- 
of -closest approach, denoted by; EL = 25° South, (25 Point of Closest approach 
to Wallops (W) and EL = 25° North (25N). The corrections are given in Table 6.3 
in this order. 


The overall scaling factor including geographic and elevation angle dependence is 
defined by the ratios 


and 


Gj(S/W) 

= R-(25S)/R-(W) 

G.(W/N) 

= R.(W)/R.(25N) 

G.(S/N) 

= R.(25S)/R.(25N) 

where 

i = o,c,4,5,7. 


For the R data, which is based on the assumption of non -geographically scaled 
c 

profile equal to the local composite, this ratio represents only the elevation angle 
scaling. This is then used as a basis to remove the first-order elevation angle 
dependence by forming the geographic scaling factox's 

F.(S/W) = G.(S/W)/G^(S/W) 

Fj(W/N) = G-(WAr;/G^(W/N) 

F.(S/N) = G.(S/N(/G^(S/N) 

tabulated in Table 6.4. 


Finally^ using the actual measurements as given by the SECOR data, (F^) the 
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TABLE 6.3 

RANGE CORRECTIONS (METERS) 
VARIOUS GEOGRAPHIC SCALINGS 


Date 

EL 

(deg) 

RO 

RC 

F7 

F8 

R5 

R7 

R8 

4 03 

25.01 S 

104.61 

83.28 

1.030 

1.090 

86.41 

85.78 

90.75 


67.80 

52.86 

46.90 

0.978 

1.000 

45.83 

45.87 

46.90 


25.03 N 

77.03 

83.35 

0.946 

0.910 

78.95 

78.85 

75.82 

4 05 

25.07 S 

65.20 

57.54 

1.112 

1.056 

64,07 

62.25 

60.76 


69.49 

36.69 

32.41 

1.013 

1.000 

32.85 

32.73 

32.41 


25.03 N 

50.77 

57.64 

0.961 

0.945 

55.56 

55.69 

54.46 

410 

24.98 S 

66.93 

64.19 

1.067 

0.944 

68.75 

68.49 

60.58 


82.53 

37.26 

33.73 

0.995 

1.000 

33.55 

33.56 

33.73 


25.03 N 

72.83 

64.24 

0.954 

1.060 

61.36 

61.29 

67.95 

4 12 

24.99 S 

54.50 

55.06 

1.115 

1.038 

61.79 

61.39 

57.19 


56.19 

32.24 

34.12 

1.011 

1.000 

34.52 

34.49 

34.12 


25.04 N 

51.80 

55.11 

0.948 

0.962 

52.37 

52.24 

53.05 

4 17 

25.07 S 

70.95 

62.03 

1.106 

1.103 

68.75 

68.60 

68.44 


81.74 

35.59 

32.90 

1.003 

1.000 

33.10 

33.09 

32.99 


25.00 N 

60.05 

62.34 

0.945 

0.899 

59.04 

58.91 

56.06 

5 24 

25.24 S 

125.04 

88.24 

1.190 

1.088 

106.57 

105.01 

96.01 


84.28 

51.44 

48.69 

0.977 

1.000 

47.45 

47.57 

48.69 


30.20 N 

42.25 

81.03 

0.920 

0.752 

74.79 

74.55 

60.93 

5 25 

30.30 S 

52.90 

47.12 

1.334 

1.016 

62.95 

62.86 

47.86 


76.49 

31.31 

28.42 

1.105 

1.000 

31.66 

31.40 

28.42 


34.88N 

48.77 

43.44 

1.103 

1.008 

48.30 

47.91 

43.79 

5 29 

25.33 S 

83.72 

62.83 

1.115 

0.975 

70.92 

70.06 

61.26 


62.18 

46.67 

37.42 

0.960 

1.000 

35.40 

35.92 

37.42 


25.22 N 

77.09 

62.23 

0.848 

0.951 

53.85 

53.62 

60.14 

5 30 

25.11 S 

85.14 

54.82 

1.170 

0.893 

64.35 

64.14 

48.95 


79.51 

37.62 

29.81 

0.948 

1.000 

28.25 

28.26 

29.81 


34.94N 

54.15 

45.61 

0.897 

0.924 

41.16 

40.91 

42.14 

6 04 

32.26 S 

62.31 

46.01 

1.125 

1.037 

52.09 

51.76 

47.71 


59.61 

39.50 

32.01 

0.938 

1.000 

30.20 

30.03 

32.01 


25.18 N 

59.10 

52.96 

0.807 

0.950 

43.15 

42.74 

50.31 

6 05 

25.38 S 

93.34 

67.26 

1.375 

1.008 

90.93 

92.48 

67.80 


75.16 

39.40 

36.88 

0.974 

1.000 

35.96 

35.92 

36.88 


35.18 N 

56.29 

55.81 

0.809 

0.916 

45.82 

45.15 

51.12 

6 11 

39.96 S 

41.18 

43.40 

1.123 

0.995 

47.98 

48.74 

43.08 


68.99 

31.35 

31.87 

1.061 

1.000 

33.60 

33.81 

31.87 


25.13N 

65.55 

57.78 

1.076 

1.050 

62.14 

62.17 

60.67 

6 13 

39.88 S 

46.49 

48.57 

1.021 

1.002 

49.60 

49.59 

48,67 


77.97 

32.74 

34.25 

0.979 

1.000 

33.57 

33.53 

34.25 


25.33 N 

76.42 

63.72 

1.021 

1.186 

65.65 

65.06 

75.57 

6 21 

30.18 S 

47.20 

33.71 

1.771 

0.993 

56.69 

59.70 

33.47 


49.68 

30.82 

24.71 

1.117 

1.000 

27.23 

27.60 

24.71 


45.13 N 

30.82 

26.25 

0.951 

0.978 

25.12 

24.96 

25.67 

6 25 

29.82 S 

65.19 

50.69 

1.657 

0.981 

81.98 

83.99 

49.73 


83.98 

33.84 

29.52 

1.020 

1.000 

29.95 

30.11 

29.52 


39.99N 

44.39 

42.46 

0.751 

0.957 

32.69 

31.89 

40.63 
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fractional error is taken as 

F.(X) - F^OO 


efiO = 




fs/w 

X =<W/N 

\sm 

The mean and estimated standard deviation of these errors are given in Table 6.4, 


Within the limitation of the small data sample here represented^ the following con- 
diisions emerge. 

1) In all cases - even unsealed j die mean scaling factor errors are relatively 
small compared to the random part. There is no significant indication 

of bias. 

2) Without any scaling at all the standard deviation of the error is of the order 
of 13% (S/W), 28% (W/N), 24% (S/N). 

3) The 5, 7 corrections which are based on the predictions are no signi- 
ficant improvement over no geographic scaling at all. 

4) The best scaling is #8 which is the linear scaling in latitude based on 

bottoms ide f F« measurement, 
o z 

This conclusion should be further qualified by noting that these passes were all 
essentially South -North and represent only two times of day. 
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TABLE 6.4 

GEOGRAPHIC SCALING FACTORS, F(X,Y) 


DATE - 1968 



U W 
Ph 

tl< Ch 


STD. 



4-03 

4-05 

4-10 

4-12 

4-17 

5-24 

5-25 

5-29 

5-30 

6-04 

6-05 

6-11 

6-13 

6-21 

6-25 

MEAN 

DEV. 








F(S/W) 










SECOR 1.114 

1.000 

0.943 

1.047 

1.055 

1.341 

1.019 

1.068 

1.231 

1.097 

1.299 

1.153 

1.255 

1.123 

1.122 

0 

0 

RC 

1.000 

1.000 

1.000 

1.000 

1.000 

l.COO 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

-0.102 

0.133 

R5 

1.062 

1.098 

1.076 

1.109 

1.105 

1.239 

1.199 

1.193 

1.239 

1.200 

1.387 

1.020 

1.051 

1.526 

1.594 

0.079 

0.173 

R7 

1.053 

1.071 

1.076 

1.102 

1.103 

1.218 

1.207 

1.162 

1.234 

1.199 

1.412 

1.014 

1.043 

1.586 

1.624 

0.079 

0.185 

R8 

1.090 

1.055 

0.943 

1.038 

1.103 

1.088 

1.016 

0.975 

0.893 

1.037 

1.008 

1.050 

1.186 

0.993 

0.981 

-0.076 

0.120 








F(W/N) 









[ 

SECOR 

1.222 

1.284 

0.974 

1.005 

1.120 

0.494 

1.019 

0.978 

0.941 

0.904 

0.944 

0.965 

1.001 

0.941 

0.912 

0 

0 

RC 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.064 

0.280 

R5 

1.031 

1.051 

1.041 

1.064 

1.059 

0.947 

0.998 

0.900 

0.952 

0.864 

0.842 

1.049 

1.042 

0.868 

0.759 

0.020 

0.255 1 

i'.7 

1.034 

1.045 

1.043 

1.066 

1.062 

0.942 

0.998 

0.883 

0.946 

0.860 

0.831 

1.059 

1.043 

0.851 

0.736 

0.014 

0.255 1 

R8 

1.099 

1.058 

0.945 

1.038 

1.112 

0.752 

1.008 

0.951 

0.924 

0.950 

0.915 

0.995 

1.002 

0.978 

0.957 

0.022 

0.147 

i 








F(S/N) 









1 

SECOR 

1.359 

1.286 

0.919 

1.053 

1.187 

2.718 

1 000 

1.093 

1.308 

1.214 

1.376 

1.196 

1.253 

1.193 

1.230 

0 

0 ] 

RC 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

-0.182 

0.238 ^ 

R5 

1.095 

1.155 

1.121 

1.181 

1.170 

1.308 

1.202 

1.325 

1.301 

1.390 

1.647 

0.973 

1.009 

1.757 

2.101 

0.071 

0.297 

R7 

1.089 

1.119 

1.118 

1.176 

1.170 

1.293 

1.210 

1.315 

1.304 

1.394 

1.700 

0.958 

1.000 

1.863 

2.206 

0.081 

0.325 

R8 

1.198 

1.117 

0.892 

1.079 

1.226 

1.337 

1.008 

1.025 

0.966 

1.092 

1.101 

1.058 

1.184 

1.015 

1.025 

-0.119 

0.172 


On 

• 

H- 
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7. MULTIPLE FREQUENCY IONOSPHERIC MEASUREMENTS COMPARISONS 

Two of tile systems used in the WICE made multiple frequency measurements 
of the ionosphere to provide their own corrections ^ but incidentally provided a 
valuable basis of intercomparison and reference for other eictemal ionospheric 
corrections. 


The SECOR system utilizes a pair of 2:1 coherently related carriers on the 
down -link at 224.5 and 449.0 MHz. Group range is measured on each and 
the difference provides an absolute measurement of the ionospheric error as 
explained in section 2.3.1. 


The TRANET system transmits from the satellite pure CW frequencies of 162^ 
324 and 972 MHz in a coherent 1:2:6 relationship. Processing of the data as 
explained in 2.4.4 yields an estimate of the ionospheric doppler or phase rate 
error Afj. 


The SECOR and TRANET measurements can be scaled to any reference fre- 
-2 

quency using the f law; the standard frequency for these comparisons has 

been the SECOR equivalent frequency, i.e. , the single frequency equivalent 

on which would be observed the same errors as are seen on the actual primary 

up and down link frequencies defined by 

"~2 _ 1 . 1 ‘ 

f2 j2 |2 

eq up down 
or 



(7.1) 


(7.2) 


(7.3) 


(7.4) 
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The SBCOR measurements are referred to f by 

eq ■' 

\2 


^secor 


eq 


and the TRANET measurements are referred to f „ by 

eq 

f, 


AR 


eq, 


^tranet 


’T nilV^eq 


where ARg and Afj are as defined in sections 2,3.1 and 2.4.1, 


The fundamental problem remains that the TRANET measurements are 
basically range-rate error whEe the SECOR measurements are in range. In 
principle the comparison could be made either in die rate domain by differ- 
entiating SECOR data or in the range domain by integrating TRANET^ but 
with the complication of an undetermined constant of integration. Both have 
been tried and it has emerged that the range comparison is more instructive, 
and that the undetermined constant can be resolved by a suitable regression 
procedure. 


(7.5) 


(7.6) 


The general approach to this regression may be described as follows ,.,nd is 
described in more detafl in Refs. 29 and 30. 


The measured ionospheric error is modelled as 


AR*(t) = AR'^CO+B + e 


(7.7) 


where e is a randomfitedng error to be minimized 

B is a bias term to be recovered 
c 

AR (t) is the computed range error at time t modelled as either 
1) "Time Gradient" model 

AR"(t> , [N„ + N„(t.^y]AB,(t) 


(7.8) 
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2) "Latitudinal Gradient" model 

AR‘=(t) = [N^ + N^(L(t)-L(y)]ARj(t) 

3 

where = maximum ionization density (e/m ) 

gradient of max. ionization density, (e/m /sec) 

/N 3 

= latitudinal gradient of max. electron density (e/m /deg) 
t^ = time of closest approach 

L = latitude 

ARj^(t) = per unit maximum density range error at time t 

^sat 

s J (n (h(s)) - 1) ds 


(7.9) 


where 


h(s) is taken as the straight path 


h(s) = + s^ + 2as sin E(t) ' - a 

a = radius of earth 

n(h) - 1 = - 40. 25 (^)/4z 

Nj^(h) is the per unit maximum density ionospheric model taken in 
various cases as either "Chapman" or "empirical" 
"Empirical" model defined by Figure 7.1 
"Chapman" model defined by 


NjO.) 

= exp(l-z-e"^) 
h-h 

<e/m^) 

(7.13) 

z 

B 

11 


(7.14) 


= 375 km 

(nominal) 

(7.15) 


= 1 [30 + 0. 2 (h^ - 200)] (km) 

(or in some experiments 

K 

(7.16) 


Hg = -^ [30 + 0. 2 (hjj^ -200)] km) 

(7.17) 
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As a first test of this regression procedure it was excercised on the problem 
of determining bias if any in the SECOR two -frequency ionospheric data 
aCOR) (Ref. 30). 


Several preliminary runs with variations in the model by adding an N term 

showed that the regression was badly unstable in this case (overmodelled). 

« 

Similarly several runs eliminating the N term showed that the gradient term 
was necessary in general to avoid badly biased resixlts. Essentially all the 
succeeding r.'i's then were with the three unknown terms as described, B, 
and N_ (or N^^). 


m 


m 


Next a series of runs were made to test the sensitivity of the results to the 
ionospheric modelling, parameters, and h„ (in eqs. 7.14 and 7.17). 

The results are shown in Table 7.1. In interpreting these results it would be 
hoped that the recovered bias, B, would be relatively insensitive to the 
modelling. In fact, it is found that over the reasonable range of variation of 
the ionospheric layer shape parameters, the recovered bias varies by as much 
as 3. 4 meters but generally less than 1/2 meter. The recovered and 
on the other hand should be expected to vary directly with the modelling para- 
meters in such a way as to tend to hold the total ionospheric content more or 
less fixed. From equations(7.l3) to (7.17) the totallayer refractivity content, 
, which is just the range error on a vertical ray is given for the 434.26 MHz 
SECOR equivalent frequency by 

o' s^ m, m 


Similarly 

M ,h„,N 1 
o' s^ m’ m^ 


= J (n--l) dh 

o 


40.3 

T 


s m 


40.; 

3e200 ^ ^ . r. 

(434.26)^10^'^ 3 ”^e/m 

3.873 • 

10"^ k [h - sol 

°^e/cc "L J 

3.873 * 

10"^ k Ih 

^e/cc/min ^ ^ ’^km 
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TABLE 7.1 

SECOR BIAS REGRESSIONS 

TESTS OF SENSITIVITY TO LAYER SHAPE PARAMETERS 


FIXED 

LAYER 

SHAPE 

PARAMETERS 


RECOVERED 

PARAMETERS 


INFERRED SECOR 

MOMENTS (RAW) 

(Should be Invariant (Adj. to 

to Model for Date) Vertical) 


Date 

k 

(scale 
ht factor) 
in eq. (7.17) 

h 

m 

Used 

SECOR 

BiaSjB 

(meters) 

N^(el/cc) 

N„ 

m 

(el/cc/min) 

M 

o 

(meters) 

M 

o 

(meters /min) 

M 

o 

(meters) 

04/10/68 

5.24 

350 

-2.3 

0.6216 X 10^ 

6574 

37.79 

.40 

36.12 

04/10/65 

5.0 

350 

-2.1 

0.6603 


6858 

38.31 

.40 

04/12/68 

5.24 

350 

-3.9 

0.4998 


- 3502 

30.39 

- .21 

27.44 

04/12/68 

5.0 

350 

-3.8 

0.5196 


- 3640 

30.15 

- .21 


04/17/68 

5.24 

350 

-2.1 

0. 5990 


-12919 

36.42 

- .79 


04/17/68 

5.0 

350 

-2.0 

0.6194 


-13450 

35.94 

- .78 

3S • 27 

04/05/68 

5.24 

350 

+3.9 

0.4892 


-12980 

29.74 

- .79 

33,49 

05/24/68 

5.24 

375 

+ 8.3 

0.6415 


-72301 

42.26 

-4.76 


05/24/68 

5.24 

350 

+ 9.7 

0.6720 


-77546 

40.86 

-4.71 

51.40 

05/24/68 

5.24 

400 

+6« 3 

0.6222 


-67813 

44.14 

-4.81 


05/25/68 

5.24 

400 

-2.0 

0.4618 


- 3135 

32.76 

- .22 


05/25/68 

5.24 

375 

-1.3 

0.4869 


- 3307 

32.07 

- .22 

31.47 

05/25/68 

5.24 

350 

-0.5 

0.5163 


- 3507 

31.39 

- .21 


05/25/68 

5.0 

350 

-0.4 

0.5368 


- 3664 

31.14 

- .21 


05/30/68 

5.24 

350 

-4.3 

0.6280 


-12539 

38.18 

- .76 


05/30/68 

5.0 

350 

-4.2 

0.6575 


-13145 

38.15 

- .76 

uO* 

04/03/68 

5.24 

350 

+2.4 

0.7813 


-26620 

47.51 

-1.62 

AQ 1 A 

04/03/68 

5.0 

350 

+ 2.7 

0.8138 ' 


-27821 

47.21 

-1.61 
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The results of these calculations are given under tlie columns "inferred 

moments". Note tliat except for the anomolous date of 5/24/68 where ex- 

* 

tremely large gradients were observed, the results of tlie and 
calculations are indeed reasonably invariant (< 5%) to the assumed layer 
shape parameter variations. Finally in the last column are listed the SECOR 
(ICOR) measurements of adjusted to the vertical as in Table 6.2. These 
can be compared to the inferred and the inferred ICOR bias. 


Table 7.2 shows the results of a corresponding adjustment of the integrated 
TRA.NET data in column @ , and the comparison to the SECOR data in 
terms of the quantity 


= AR 


secor 


tranet 

secor 


X AR 


tranet 


in column (6) . This can be compared to the SECOR bias adjustment defined 
from (7. 7) as 

B = -AR'^ - e 

secoi' 

where AR^ is the computed, modelled range error. 


Five cases in the two sets of data overlapped and these are listed in column . 
While there may be some hint of a constant bias difference of several meters 
between and B the small data sample precludes ascribing too much 
significance to it. 


Column (E; lists the standard deviation of tlie difference between AR(t) 

^ 9 ' '^secor 

and ARj.j^j^g^ x -p over the duration of the common track after removing 

secor . 

the mean difference. The differences run in the area of 1/2 to 1 meter and when 
plotted appear to show no consistent systematic trends. 


The nature of these differences is further indicated in Figures 7.2 - 7.11 which 
show the variation over the observed pass of the difference between the SECOR 
and the TRANET observed ionospheric range error, referred to SECOR equi- 
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TABLE 7.2 


LEAST -SQUARES SOLUTION USING TRANET DATA 



A 

B 

C 

D 

f— ■■ 

, E 


SECOR*^^ 

TRANET'^* 

SECOR^-^^ 

TRANET 

SECOR^^^ 

I 

< 

®4 

DATE 

AR(t„) 

AR(to> 

^4<'o> 

BIAS 

B 

MA) 

Meter 

Meter 

Meter 

Meter 

, Meter 


4/05 

35.1 

28.4 

6.7 

3.9 1 

.60 

4/10 

35.7 

35.6 

.1 

-2.3 

.61 

4/11 

24.7 

28.6 

-3.9 


1 .55 

4/12 

32.5 

32.8 

- .3 

-3.9 

.43 

4/16 

45.6 

43.5 

2.1 


.43 

4/17 

35.0 

33.4 

1.6 

-2.1 

i 1.02 

4/18 

29.5 

31.4 

-1.9 


.43 

4/21 

40.0 

39.2 

.8 


.50 

4/23 

42.3 

34.4 

7.9 


.85 

4/26 

63.7 

68.0 

-4.3 


.70 

4/29 

34.6 

31.1 

3.5 


1 .40 

4/30 

29.4 

27.7 

1.7 

j 

i .53 

5/01 

29.0 

26.8 

2.2 


1 .40 

5/02 

44.1 

47.9 

-3.8 


.59 

5/03 

37.2 

41.9 

-4.7 


! .36 

5/04 

36.2 

36.3 

- .1 


1 .37 

5/07 

44.8 

45.1 

- .3 


1 .36 

5/08 

20.0 

27.2 

-7.2 


1 1.12 

5/25 

31.0 

26.7 

4.3 

= 0.232 m 
CT = 3.93 m 

• 

-0.5 

' .93 

1 

r 

1 

i 

! 


(1) TRANET Data is self adjusted for Bias. N^, N^ and scaled to f 

•’ ’ m' m secor 

(2) SECOR Unadjusted, at f 

S0COX* 

<®>‘^4 = ^R<'o>secor-'^R<to)tranet' ^ adjusted to 

(4) From Table 7.1, using k = 5.24, h = 350 km 

' ° s ' m 

(5) t^ = time of closest approach 
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valent frequency and adjusted for the mean difference. For tlie purpose of 
this comparison it was found necessary to smooth tlie SECOR Ionospheric data 
rather heavily by two successive passes through a 36 -second span filter, with 
rejection of all points more than 1.5o off the curve, between passes. There 
do not appear to be any consistent systematic trends in these differences . 

It is of interest to compare the results of these SBCOR-TRANET Comparisons 
with the SECOR vs Profile -Ray -Trace r'esults in Table 6.2a), which indicate 
comparable RMS differences of the order of 3 -4m. It is significant then 
that the TRANBT data can be integrated with the integration constant resolved 
internally by absolute ionospheric error estimates to at least the same order 
of accuracy as the best profile determinations . 
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FIGURE 7.2 

VARIATIONS IN THE DIFFERENCE 
BETWEEN SECOR AND INTEGRATED 
TRANET IONOSPHERIC RANGE 
ERROR MEASURBMFNT55 
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b) 


FIGURE 7.3 

DIFFERENCE IN VARIATIONS AS 
OBSERVED IN SECOR DATA AND 
INTEGRATED TRANET DATA 
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FIGURE 7.4 

DIFFERENCE IN VARIATIONS AS 
OBSERVED IN SECOR DATA AND INTEGRATED 
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FIGURE 7.5 


DIFFERENCE IN VARIATIONS AS 
OBSERVED IN SECOR DATA AND INTEGRATED 
TRANET DATA 
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FIGURE 7.6 

DIFFERENCE IN VARIATIONS AS 
OBSERVED IN SBCOR DATA AND INTEGRATED 
TRANET DATA 
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FIGURE 7.7 

DIFFERENCE IN VARIATIONS AS 
OBSERVED IN SECOR DATA AND 
INTEGRATED TRANET DATA 
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FIGURE 7.9 

DIFFERENCE IN VARIATIONS AS 
OBSERVED IN SBCOR DATA AND 
INTEGRATED TRANET DATA 
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FIGURE 7.10 

DIFFERENCE IN VARIATIONS AS 
OBSERVED IN SECOR DATA AND 
INTEGRATED TRANET DATA 
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FIGURE 7.11 

DIFFERENCE IN VARIATIONS AS 
OBSERVED IN SECOR DATA AND 
INTEGRATED TRANET DATA 
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8. SECOR-ORBIT COMPARISONS 

Reference orbits were available for the 17 ionospheric test case passes. 

These reference orbits were derived from short arc orbital fits to observed 
Laser, C-Band Radar, and optics data. Over the Wallops observation span 
the orbital data is believed to be accurate to roughly + 1 meter in the range 
component and + 1 arcsecond in the angular components. 

The SECOR Data, corrected for Tropospheric errors as described in section 2 
and for Ionospheric errors by its own two -frequency determination was com- 
pared to the orbital data. The residual data was further adjusted for SECOR 
range bias and time bias to best fit the orbital data. The adjusted residuals 
are plotted in Figures 8.1 - 8.15, and the adjustment noted. 

On the same plots for each day and to the same scale are plotted for comparison 
the ionospheric error corrections as determined by 

a) SECOR itself , 2 frequency unadjusted 

b) Raytracing of the local (Wallops) WICE composite ionospheric profile 

c) Raytracing of the geographically-scaled composite profile, using the type 

8 gradient scaling discussed in section 6.6, i.e. , scaling proportional to 
^0^2 latitudes of Grand Bahama, Wallops, and Ottowa, 

with linear interpolation for intermediate latitudes. 

Comparing the geographically scaled versus the unsealed composite raytrace 
computations to the SECOR ionospheric measurement confirms the general con- 
clusion of section 6.6, namely that the best of the predicted latitudinal gradient 
modifications is not very good, in fact, hardly worth while. 

Comparing the measured ionospheric correction with the orbital residuals after 
correction, plotted to the same scale, indicates the high degree of correction 
achieved for the SECOR data. 

The recovered biases in range and time for the various passes are listed 
together in Table 8.1. While the overall bias trends are clearly persistent and 
statistically significant, it cannot necessarily be inferred that the discrepancies 
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4; 


noted are in fact correctly attributed to SECOR range error or timing error. 
The recovered timing biases are in fact one or two orders of magnitude 
larger than should be considered reasonable and can probably be ascribed to 
aliasing of second-order (i.e., non-linear) latitudinal variations into recovered 
time bias. 
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TABLE 8.1 

BIAS ADJUSTMENTS; SECOR-ORBIT FITS 


Date 

Recovered 

Bias 

Adjustments 

Range Time 

Bias Bias 

(MiUi- 

(Meters) Seconds) 

RMS 

Adjusted 
Fit Error 

(Meters) 

4/03 

-11.12 

-1.06 

1.85 

4/05 

-16.15 

-1.16 

2.41 

4/10 

-10.19 

-0.82 

2.05 

4/12 

-14.30 

-0.41 

2.00 

4/17 

-13.21 

-0.92 

1.81 

5/24 

-20.8 

-2.31 

1.81 

5/25 

-13.5 

-1.76 

1.79 

5/29 

5/30 

-15.1 

-1.16 

1.67 

6/04 

-19.4 

-0.87 

1.50 

6/05 

-17.7 

-0.44 

1.85 

6/11 

-19.4 

-0.87 

2.16 

6/13 

-19.1 

-0.96 

1.58 

6/21 

-22.7 

-1.39 


6/25 

-10.1 

-1.37 



Mean -15.91 -1.11 


I RMS 4.05 .50 

I 




Ionospheric ranoe error (meters) 
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9. SUMMARY CONCLUSIONS 


Ray Trace (Section 3) 


S.I) After minor modification of the density extrapolation routine above the 
top of the ionosphere the REEK Program was found to be an extremely 
and consistently accurate basis for ray trace comparisons. 


3.2) For certain purposes it is a reasonable approximation to ignore bending 
in the path integrals. For a simple slab model ionosphere the error in 
so doing is approximately 


2 

(1 - 1) tan^5 sec? 

where N = (average) refract ivity 

T = layer thickness 
h = satellite height 
? = zenith angle 

The error is quadratic in N and vanishes at vertical incidence. 


3.3) For the same slab model it is shown that the group range error is given 
exactly by the integral of group refractivity along the phase path. 

3.4) Numerical comparisons of REEK and straight line raytrace on actual 
layers confirm the above general relationship. 


3.5) For a representative ionospheric case of MHz, 

f = 434 MHz (SECOR), « 84 • 10 the maximum error in the 

straight line assumption is 0.11 meters out of 73 meters. 


3.6) For the troposphere, the range error in ignoring bending becomes 

significant, i.e. , equal to or greater than 2% for elevation angles less 
than about 10°. This is a significant limitation of any troposphere 
correction method which ignores bending. 


( > 
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_3 

3.7) The neglect of f (gyromagnetic) terms in the ionospheric refractivity 

can result in errors as great as 2% at frequencies of the order of 100 MHz. 

As this can be compared to the neglect of bending error, which varies as 
-4 

f , and is typically 0.2% at 100 MHz. Thus above 100 MHz the neglect 
of bending error is generally less than 1/10 that due to neglect of f ternis 
in tlie refractivity. 

3. 8) It is proven analytically in Section 3. 3 that the ionospheric group range 
error is given exactly by the integral of group refractivity along the phase 
path. 

3. 9) The error in the use of the superposition principal to determine ionospheric 
and tropospheric errors separately was studied numerically. With typical 

_3 

values at 100 to 400 MHz the error is less than 10 of the range error at 
low elevation angles. 


Analytic Corrections (Section 4) 

4.1) All of the tropospheric analytic forms considered are of essentially 
identical form above 10° elevation. 

4.2) Below 10° elevation angle the NAP-1 and GDAP formulations are a good 
approximation to the correct variation down to the lowest angles . These 
corrections are of the general form 


AR(E) = AR(90°)cscE 

2 


1 + a/i + .004 csc^E 

AE(E) - ctn E r 

' 1 

s L 

/ — z": — T~' -I 

1 + A 

/l + . 004 CSC E 


4. 3) A special regression study was conducted, based on 85 days of actual 
radiosonde profile data taken during WICE and considering various 
possible models for the total range error in terms of measured values 


I 

i 

I 

I 

i 

1 


! 
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of pressure, temperature, and humidity at ground level. The best model 
was found to be the one which depended on Ni^P, T, H), i.e. , the ground 
level refractivity alone. This is the same in form as the NBS regression 
model. 

4.4) The predictive capability of this special regression for the 85 day data 

from which it was derived was compared to that of the t vo NBS regression 
models established on Cape Canaveral data and on a widespread US data 
base, and to the Hopfield bi-quartic prediction. The standard of com- 
parison was REEK raytracing of the radiosonde measured profiles. 

The special Wallops regression was of course best (.149 m rms) followed 
by Hopfield (.223 m rms) NBS "Cape Canaveral" regression (.258 m rms), 
and NBS "Standard Sample regression" (.281 m rms). 

' 4.5) For the Ionosphere the minimum elevation angles are never less than 
about 18° in the ionosphere even for horizontal takeoff at the ground. 

For this reason, most of the analytic forms studies, with the exception 
of the Freeman and GEOVAP formulations provide reasonably good 
models down to the lowest ground elevation angles. 


Moment Expansions (Section 5) 


5.1) A moment series expansion, useful for ionospheric and tropospheric range, 
elevation angle, and range rate errors is devel oped in the form 


M 

AR(E) I G^(E)M^ 
m=0 

til 

where = m moment of the profile refractivity distribution 

= fN(h)(h-h^)"'dh 
o 


- 9.3 - 


Communications Research Laboratories 




I 


GJE) = 

-- "c 

The utility of the expansion derives from the fact that the layer charac- 
teristics are totally characterized by the moments (independent of E) 
while the elevation angle dependence is totally characterized by the coef- 
ficients Gj^(E) (independent of the layer characteristics). 

.. - . . i » 

5.2) For most purposes the expansion need be carried only to the M = 2 
term for full accuracy. 

5. 3) The series is always convergent for the ionosphere but is convergently 
beginning but ultimately diverges for the troposphere below about 2 
elevation angle. 

5.4) For both the ionosphere and the troposphere the moment series expansion 
corrections for AR, AB, and AR agreed better with the REEK ray- 
traces than did any of the other analytic corrections tested. 

Profile Determinations (Section 6) 

6.1) Four different methods of profile prediction were available for comparison, 
the Jackson composite bottoms ide -topside sounder reductions, the ESSA 
bottoms ide sounder extrapolations, the Freeman model based on 3 -month 
predictions (no longer published) and based on long-time predictions 
(CCIR Report 340) of f^F 2 and f 3000 MUF. 

6.2) It was found that the discrepancies between the various profile estimates 
could be expressed efficiently in terms of the moment series particularly 
since only the zeroeth order moment differences are significant. The 
first and second order moment differences are responsible for less than 
1/10 of the range error differences (in meters) of that due to the zeroeth 
order moment even for the worst case of lowest elevation angles 
(Figure 6.2). 

6.3) The four predictions listed above were compared in terms of ' th 

the SECOR measurements as a reference (Table 6.2). The rms differences 


geometrical coefficient 

ih sec (5(h)) 

dh , _ 
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relative to SECOR and therefore rough estimates of the absolute error of 
the various predictions were 


Profile Measurements 

m 

a 

rms 

meters 


Jackson Composite 

-4.1 

2.5 

4.8 

13 

ESSA Bottomside 
Extrapolation 
(only 6 cases) 

-1.5 

3.4 

3.7 

10 

Profile Predictions 





CCIR Report 340 

-5.2 

7.1 

8.8 

24 

3 Mos. Advance 
Predictions 

-9.1 

6.9 

11.4 

32 


6.4) The SECOR measurements demonstrated significant north-south gradient 
effects. Several methods were studied to estimate these geographical 
gradients in terms of both measured and predicted profiles. From 
measured profiles, gradients were inferred from the measured dif- 
ferences between Ottowa, Wallops andGrand Bahama. From the pre- 
dicted profiles, gradients were inferred from the differences of predicted 
densities at each point along the ray, as a function of the actual geograph- 
ical coordinates of the point on the ray. 

None of the geographical scalings were particularly effective in predicting 
the actual (SECOR) north-south gradient effects. The most effective was 
that based on linear scaling in latitude based on bottomside f F, measure 
ments from Ottowa, Wallops and Grand Bahama. For the 17 WICE ionos- 
phere sample days tlie average south/north assymetry factor was 1.29,*. 1. 
Application of the bottomside f^F 2 scaling predicted the assymetry with 
an rms residual of .17 which represents a 60% unpredictable residual. 


Multiple Frequency Comparisons (Section 7) 

7.1) The SECOR and TRANET systems botli made internal measurements of 
ionospheric error, the former in terms of range and the latter in terms 
of range -rate. In order to compare the two, either SECOR data had to be 
differentiated or TRANET had to be integrated, with the resulting constant 
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of integration determined by other means. Regression methods were 
developed for the estimation of that constant of integration and applied to 
both SECOR and TRANET. 

7.2) The first series of tests of the regression procedure were in an attempt 
to determine whether there was any systematic bias in the SECOR 
ionospheric measurements . 

The ionosphere was modelled as of a given shape (Chapman or other) 
with fixed thickness parameter and undetermined maximum density at 
point of closest approach and undetermined north -south (or time) 

gradient (Njj^). 

The recovered biases (Table 7.1) ranged from about -4 to +9 meters 
(at 434 MHz). On any particular day the biases were reasonably insensitive 
to modelling assumptions (thickness and height of maximum) but from day 
to day there was no significant persistence of bias. It can probably be 
concluded that the recovered biases do not represent real SECOR errors 
but regression modelling errors due to the fact that the ionosphere de- 
parted significantly from a simple linear north -south gradient. 

7.3) The TRANET data were integrated and a similar regression performed 
to fix the constant of integration. The results were compared to SECOR 
in terms of range error at point of closest approach (e^). The mean dif- 
ference (Table 7.2) was 0.232 m with a standard deviation of 3. 93 m or about 
10% (at 434 MHz) which is closely comparable to the differences Detween 
SECOR and the best profile determination ray trace results. This provides 

a valid measure of accuracy with which absolute ionospheric measurements 
can be recovered from integrated TRANET by regression for the constant 
of integration. Point by point comparison of SECOR vs integrated TRANET 
yielded a standard deviation of 0. 59 m which is a measure of the random 
instrumental disagreement of the two systems without regard to the con- 
stant of integraticii. 
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SECOR-ORDIT Comparisons (Section 8) 

8.1) SECOR range data was compared to the best short arc orbital data as 
defined by LASER, Camera, and C-Band radar. Also in the same com- 
parisons, the SECOR ionospheric data was compared to ray tracing 
ionospheric corrections based on the Jackson composite data and based 
on the geographically scaled Jackson composite data utilizing the scaling 
method previously found best, i.e. , linear scaling in latitude based on 
bottomside f^F 2 measurements at Ottowa, Wallops and Grand Bahama. 

8.2) The ionospheric data comparisons reiterate the same conclusion found 
previously that the best geographic scaling is not very good at accounting 
for actual observed variations of ionosphere along the trajectory since 
the latter do not usually fit a constant gradient model very well. 

8.3) The SECOR vs ORBIT range data comparisons yielded a statistically 
significant mean bias of -15.9 meters with a day-to-day standard 
deviation of 4.05 meters and a much larger than expected time bias 
of -1.11 milliseconds with o = .5 milliseconds. These biases are 
probably largely a result of aliasing of non-linear variations of actual 
ionospheric content into the recovered time bias term. 
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